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Abstract 


The simulation of ultrasonic wave propagation is one of the emerging fields of research in 
non- destructive testing & evaluation Looking at the underlying physics, the challenging 
problems in ultrasonic NDT are on one hand, caused by micro and macrostructural material 
characteristics as inhomogeneity and anisotropy, and on the other hand -by the presence of 
complex component defect geometries. Both can lead to considerable problems in testing 
performance and in the interpretation of inspection results. Straightforward solutions like 
experience-based adaptation in the experimental set-up are in general insufficient, insecure 
and cost-intensive. The better choice is to combine practical experience with theoretical 
modeling. 

In ultrasonic nondestructive testing (NDT) use is made of the properties of elastic waves 
propagation in solids in order to detect the defects and materials inhomogeneities. However, 
modem materials mostly exhibit anisotropic elastic behavior leading to complicated wave 
propagation phenomenon. To ensure the reliability of ultrasonic inspection techniques, these 
material properties as well as the influence of microstractural inhomogeneities and the effects 
of non- planar surfaces and interfaces on ultrasonic wave propagation have to be taken in to 
account. 

In the present work a numerical model has been developed for understanding the ultrasonic 
wave propagation in composite materials. Typically wave speed, energy flux deviation have 
been the primary parameters. From this model direction of low inspection sensitivity and also 
the regions of material to which no ultrasound penetrates can be predicted The relevance of 
this work to the ultrasonic non-destructive testing of composite material has been discussed, 
concluding that care is needed over the choice of wave modes, probe angels and receiver 
position, to ensure sensitive inspection of the whole component 
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Chapter 1 

Introduction 


1.1 Introduction 

The search for ever lighter, stiffer, and tougher substances from which to construct aerospace 
vehicles has led inexorably towards more highly engineered, lower density materials like 
composite materials. The novelty and expense of composite materials has meant, historically, 
that they are thoroughly inspected prior to use, typically using ultrasonic non-destructive 
testing. 

Ultrasonic materials characterization is the most important application of ultrasonics in 
aerospace engineering and engineering mechanics. Historically, ultrasonic nondestructive 
testing (NDT) has been used almost exclusively for detecting macroscopic discontinuities in 
structures after they have been in service for some time. It has become increasingly evident 
that it is practical and cost effective to expand the role of ultrasonic NDT testing to include 
all aspects of materials production and application. Research efforts are being directed at 
developing and perfecting NDT capable of monitoring (i) material production processes, (ii) 
material integrity following transport, storage and fabrication, and (iii) the amount and rate of 
degradation during service. In addition, efforts are underway to develop techniques capable 
of quantitative discontinuity sizing, permitting determination of material response using 
fracture mechanics analysis, as well as techniques for quantitative materials characterization 
to replace the qualitative techniques used in the past. Ultrasonic techniques play a prominent 
role in these developments because they afford useful and versatile methods for evaluating 
microstructures, associated mechanical properties, as well as detecting microscopic and 
macroscopic discontinuities in solid materials. 
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But, the strong, almost pathological, anisotropy of a composite laminate has significant 
implications for its elastic behavior and for the propagation of ultrasound. For this reason the 
use of ultrasonic waves to characterize the mechanical condition of composites has presented 
engineers inspecting these advanced materials with special challenges. The main objective in 
this thesis is to study the behavior of ultrasonic wave propagation in composite materials and 
its use in their nondestructive characterization. 

From the beginning, Ultrasonics has been the tool of choice to inspect composites, since the 
likely defects and important material properties are most easily, and inexpensively, 
uncovered in ultrasonic nondestructive testing. In ultrasonic NDE, real materials with more 
complex elastic properties (anisotropy, inhomogeneity, nonlinearity, attenuation, dispersion, 
temperature- dependence, etc.) are considered. The primary purpose of ultrasonic NDE is to 
understand the wave-material interaction and assess the desirable material properties from 
the observed deviation in the ultrasonic response from that of an ideal, defect- free medium. 

The growing need for their proper testing of composites/composite laminates raised 
considerable interest in studying elastic wave propagation in anisotropic media. The main 
problem of ultrasonic testing methods used in the field of quantitative nondestructive 
evaluation arises from the difference of wave propagation direction and direction of energy 
flow. Complex microstructure, macrostructure, or both present in many of today’s structural 
component lead to complicated ultrasound paths due to beam profile distortion and beam 
path bending. Further in composite materials the direction of wave vector and direction of 
energy flow is not the same, which lead to beam skewing effects. In this respect 
mathematical modeling has evolved as an important tool providing assisting analysis and has 
optimized experimental setups. A particular interest exists in the propagation of pulses as 
applied in ultrasonic testing, where proper understanding and modeling promises optimized 
and effective application. 

1.2 Literature Survey 

Theoretical studies on the elastic wave propagation found in literature are manifold, most of 
which dealing with isotropic media, only a few are concerned with wave propagation in 
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anisotropic material. Extending the scope of interest to the wave scattering problem in 
anisotropic materials for purpose of location and classification of defects and other 
inhomogeneities, respective studies are approximate in character and mostly restricted to 
particular geometries and materials symmetries. The theory presented is obtained by applying 
a coordinate free approach, which had been first used by Chen [1], making it useful in 
dealing with lay structures. 

Synge [2] developed the theory of slowness surfaces after that Musgrave [3] first introduced 
the idea of superimposing the slowness surface on the geometric space, which contains the 
incident and emerging waves. 

Previously, modeling techniques for inhomogeneous materials have been developed with 
emphasis on inspection of austenitic welds and claddings. Although the numerical based 
methods such as finite difference, finite element or finite integration provide good results, 
they suffer from the requirement of large computation times and hardware requirements. 
Computationally advantageous analytical ray tracing method developed in this work is very 
useful in explaining certain phenomenon and in optimizing ultrasonic inspection techniques. 

In the literature, only a few ray-tracing algorithm for wave propagation in anisotropic 
inhomogeneous materials are discussed. Silk [4] divides the austenitic weld region into a 
number of quadrilateral regions, each with a fixed direction of grain crystals. Ogilvy [5] 
defined the inhomogeneity of the weld by empirical analytic function that describes the 
orientation of grain crystals continuously. 

Leander [6] developed the relation between the wavefront speed, which is the speed of 
propagation of a front of a pulse and group velocity concept. Crandell’s [7] work was based 
on use of slowness diagrams to represent wave reflections. 

Very few quantitative results have been published on the reflection and refraction of plane 
elastic waves from an interface between anisotropic media. Musgrave and Fedorov [8] 
summarized the general theoretical results for the reflection- refraction problem in their 
comprehensive books. Auld [9] discussed several examples for planes of symmetry. Henneke 
[10] and McNiven et al [11] reviewed Fedorov's method and discussed critical- angle 
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phenomena for the elastic waves at an interface. These works mainly focused on finding 
reflection-refraction angles using geometrical plotting on slowness surfaces. Vandenbossche 
et al. [12] studied the mode transition in unidirectional composites. 

More recently, Rokhlin et al. [13] have described a unified approach to the numerical 
solution of the reflection-refraction problem for generally anisotropic media. Ogilvy [14] also 
studied the influence of austenitic weld geometry and manufacture on ultrasonic inspection 
of weld joints. Ogilvy [15,16] was the first to purpose a ray- tracing model for ultrasonic 
NDT in austenite welds. Schmitz et al. [17] modified the Ogilvy’s model and proposed a new 
model for austenite materials to follow the longitudinal, horizontal and vertical polarized 
shear wave propagation from the base material through the cladding in three dimension. 

Marker et al. [18] presented full numerical solution of wave equation and showed that ray 
tracing method gives similar and consistent results. Spies [19] studied elastic wave 
propagation in general anisotropic media using Green’s function and elastodynamic 
holography. 

1.3 Present Work 

The nature of wave propagation in anisotropic and inhomogeneous media is significantly 
different from that in isotropic media. The basic difference may be formulated in the 
following manner: 

(1) For an arbitrarily selected direction in an anisotropic material, the propagation of 
three elastic wave modes is possible. These three are pure transverse, quasi- 
transverse and quasi- longitudinal. For special directions called acoustic axes the 
velocity of two transverse waves coincide. 

(2) The displacement direction for each of the wave is uniquely determined by the 
direction of wave propagation. Only for the propagation along acoustic axes can 
the transverse wave be arbitrarily polarized as in an isotropic material. 
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(3) Each of these waves has different phase and group (ray or energy) velocities. The 
ray velocity is greater than or equal to the phase velocity and its direction does not 
coincide with wave normal. 

To ensure the reliability of ultrasonic inspection techniques, these complex material 
properties as well as the influence of microstructural inhomogeneities and the effects of non- 
planar surfaces and interfaces on ultrasonic wave propagation have to be taken in to account. 
In this respect, simulation and optimization in ultrasonic testing have gained a considerable 
importance, where mathematical modeling provides a simple method of assisting analysis. In 
this presentation a ray- tracing model is developed to follow quasi- longitudinal and quasi- 
shear wave propagation in composite materials for various phase angle. 

This dissertation presents a general study of elastic wave propagation in inhomogeneous and 
anisotropic media. For the present work transversely isotropic symmetry assumption is used, 
which is the macroscopic symmetry of unidirectional fiber composites, extruded metal matrix 
composites, as well as ideally structured columnar- grained steels. The results are most 
general in that they are derived for arbitrary orientation of the rotational symmetry axis 
which-for the sake of simplicity-will in the following to be referred as the fiber direction. The 
theory presented is obtained by applying a coordinate free approach making it useful in 
dealing with lay structures. Subsequent chapters will deal with the plane wave solution of the 
equation of motion, which is characterized by the respective wave velocity, polarization, and 
group velocity vectors, the fiber direction being a variable. Then reflection and refractions of 
plane wave at an interface between two arbitrarily oriented homogeneous transversely 
isotropic layers are considered. Finally an algorithm to numerically simulate file path of the 
refracted wave is developed. 

An analytical ray tracing method has been developed which is very useful in explaining 
certain phenomenon and in optimizing ultrasonic inspection techniques. This ray tracing 
algorithm can be used to trace the ultrasonic wave propagation from the search unit through 
the composite base material and various inhomogeneities. The algorithm calculates the 
direction of the sound field, but does not render any amplitude information. Both travel path 
and travel time information can be inferred from these model. This model can be used to 
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calculate; the slight deviation which occurs in the direction of wave vector and direction 
energy flow, and to determine the ray behavior at an interface between two regions of 
different material properties. Presently a two dimensional analysis has been done and this can 
be extended to three dimensions with some modifications. Within this region, several 
defective regions may exist, those being circular, rectangular and elliptical etc. Interfaces 
between different regions are treated as locally planar. 


1.4 Thesis Organization 

Chapter II gives the fundamentals of elastic wave propagation in anisotropic media. It 
describes the equation of motion and plane wave solution of this equation. It discusses the 
fundamental equation of wave propagation in inhomogeneous and anisotropic media called 
as Christoffel’s equation. Various terms associated with the elastic wave propagation in 
anisotropic media like phase and group velocity, slowness, energy flux and critical angle are 
discussed in this chapter. Beam divergence and Beam Skewing associated with elastic waves 
are also discussed. 

In Chapter III refraction and reflection properties of elastic waves are discussed. Snell’s law 
for isotropic has been discussed and it is extended for the anisotropic case. The use of 
slowness diagrams for determining the reflected and refracted wave properties has been 
discussed in this chapter. 

Chapter IV provides the mathematical formulation for the elastic wave propagation in 
arbitrarily oriented transversely isotropic media. 

Chapter V presents the ray-tracing model, which has been proposed to simulate the ultrasonic 
wave paths in an inhomogeneous and anisotropic media. 

In chapter VI various results of the ray-tracing model developed in the work have been 
presented and these results are compared with the experimentally obtained results. 

Chapter VII presents the conclusion and scope of future work. 
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Chapter 2 

Basics of Ultrasonic wave propagation in anisotropic 

media 


2.1 Introduction 

Ultrasonics involves the propagation of acoustic waves. Therefore, it is necessary to 
understand the basic features of propagating waves and some of the mathematical equations 
governing wave propagation. Acoustics is the study of time-varying deformations, or 
vibrations in elastic media It is concerned with material particles that are small but yet 
contain many atoms. Within each particle the atoms move in unison. Therefore, acoustics 
deals with macroscopic phenomena and is formulated as if matter were a continuum. 

Propagation of ultrasonic waves in isotropic media has been well studied and can be 
predicted theoretically with good results. The problem is defined by two constants for the 
medium, such as X and fi, the Lame constants and, the incident wave speed and direction. 
For anisotropic media the problem becomes considerably more complicated, as the elastic 
constants determining the wave propagation are then directionally dependent. 

In general, for a given propagation direction, the wave equation has three solutions with 
differing speeds corresponding to three different polarization. A plane wave solution to the 
wave equation describes the possibility of three different wave modes with different phase 
velocities. These wave modes are not in general, purely longitudinal or transverse, but quasi 
wave modes. 

The effect of anisotropy on ultrasonic wave propagation can be very marked as wave 
velocities become directionally dependent; group and phase velocities are no longer 
necessarily parallel or equal in magnitude. These effects lead to beam distortion with in the 
material, unexpected beam propagation directions, ‘favoured’ propagation directions, 



alteration in spatial profiles and, for inhomogeneous media, beam paths which become 
curved as the beam propagates through the material. It is therefore important to be able to 
predict such effects to ensure that beam energy (which travels in the group velocity 
direction) is directed to the required region for inspection purpose. 

2.2 Equations of Motion 

The static equilibrium of an arbitrary, three-dimensional solid body occupying volume V 
that has both surface and body forces requires that the resultant force vector vanish, this 
gives 


dxj 


pb=0 


( 2 . 1 ) 


The term h is the body force density vector and the result is sometimes referred to as Euler's 
first relation. This is a static result, however, and to describe the dynamical motion of a 
continuum, an inertial term accounting for the particle acceleration must be included, 


3a d’u, 

— + pb, = p—- 


dx, 


dt 2 ’ 


( 2 . 2 ) 


It is also called Cauchy's equation. If body force bj vanishes for free, unforced motion, it 
gives 


d Xj p dt 2 ’ 


(23) 


and this gives dynamical equation of equilibrium for the elastic continuum. Using 
constitutive relations Eq. 2.3 can be written as 


where Cpi 
that 


r dga . 

iJ “dx. 


= P 


dt 2 


(2.4) 


is the elastic coefficient tensor. It is symmetric and positive definite in the sense 



^ 0 

for all symmetric Hy where equality is satisfied only when ny = 0 

By inserting for the strain £ki its linearized equivalent in terms of the displacement gradient, 


Eq. (2.4) can be written as 


This can be written as 



du, 

dx k 



d du. 

L + 

dXj [ dx l 


du / 

dx k 




d\_ 

dx,dx. 



(2.5) 


( 2 . 6 ) 


(2.7) 


Equation (2.7) is the equation of motion for an infinite linear elastic medium in three 
dimensions. 


2.3 Plane Wave Propagation in Bulk Materials 

The plane harmonic wave solution of Eq. (2.7) can be expressed either in index notation or 
in vector notation as 


u(r,f) = ^pexpi(k.r-G)f) 

(2.8) 

u M m ’ t ) = AP J ^Vi{k m x m -0X) 

(2.9) 


where Ao is the amplitude of the wave, p (pj) is the polarization unit vector and to is the 
angular frequency. The wavevector k or (km) is given by 

k = kn, k = — (2.10) 

A 


where n is a unit vector in the direction of wave propagation, and X is the wavelength. The 
phase of the wave $ is given by 



(p = k-r-cot 


( 2 . 11 ) 


For points of constant phase on the wave is 

k • r - at = const. (2.12) 

Taking differential and rearran ging gives 

dr 

k — =0) (2.13) 

dt 

Substituting k from Eq. (2.10) and let J; be the projection of the wavefront position vector r 
on the direction of wave propagation (£;= n . r), as shown in Fig. 2.1, 



and the propagation of a point of constant phase proceeds with phase velocity V given by 


V = ^-= — 

dl k 


( 2 . 15 ) 


in vectorial form 



(2.16) 


V = Vn = —n. 

k 


Substituting the plane harmonic wave solution (2.9) into the equation of motion (2.7) gives 

P<A=C„, .kjk,u„ (2.17) 

Or 

c v , m n j n F m -P v2 P, = 0 »- ( 2 - lg ) 

By using p. — 8 im p m , the above equation can be rewritten as 

{C„, m nn,-pV I Sjp. = 0. (2.19) 

where 8 m is the Kroneker delta. In order to have a nontrivial solution for pfo, the determinant 
of the coefficients for the linear equation Eq. (2.18) must vanish, leading to the following 
well known eigenvalue equation, also known as Christoffel's equation, 


C ,ji m n j n ,-P vl 8 i m =°- 


( 2 . 20 ) 


if A.ii = Cijbni^n and r\ = p V 2 
Eq. (2.19) can also be written as 


(Xu - TjSnn) p m =0 (2.21) 

From the properties of Cjjim, it follows that Xu is also symmetric and positive definite. 
Therefore, all the eigen values of Xu are real and positive and their corresponding 
eigenvectors are orthogonal. 



This can be interpreted as that a for a given direction of wave propagation n there will be 
three phase velocities, and three corresponding displacement vectors will be orthogonal. 
These three waves are pure shear, quasi-longitudinal and quasi shear wave. 

Pure longitudinal mode (the polarization vector p is parallel to the propagation direction m, 
1, n are the direction cosines of propagation direction) is obtained only when p l —n i (i = 
1,2,3). Pure shear mode (the polarization vector p is normal to the propagation direction n) 
are obtained only when p l n l +p 2 n 2 +p 3 n i =0 (i.e. p n = 0) In all other cases, either quasi- 
longitudinal or quasi-shear modes can only propagate. It can also be noted that any general 
wave propagation can be split in three modes: pure shear (S), quasi- longitudinal (QL) and 
quasi-shear (QS). In a elaborate form Eq.( 2.21) can be written as 


1 

>» 

1 

to 


K, 

'Px 


'O' 

K 

K - pv 2 

^23 

Pi 

> = < 

0> 



X-pv\ 

J 3. 


0 

k. / 


where X in terms of elastic constants can be given as 

K=t 2 C uu +m 2 C l2n WC uu +2mnC m2 +2n£C ul3 +2m£C un 
X 22 =£ 2 C l2u +m 2 C,,,, +n 2 C 2323 +2mnC 2223 + 2n£C 23l2 +2 m£C 22l2 

^3=^ C I3I3 +m2 Q323 + w2C 3333 +2WwC 3 3 23 +2 ^ C 3313 + 2m^C 2 3i3 (2.23) 

\ 2 = ^ 2 C ul2 +m 2 C 2212 +n 2 C m3 +mn(C 23n +C 2213 )+n£(C U23 +C m2 )+mi(C m2 +C m2 ) 

K = ^ C m3 +m2C 2in +n2 Qm +mn ( C 2m +C m2) +n ^( C im +C m0 +n ^( C n23 +C m2) 

X 22 —i C in2 +m C 2223 +n C 3323 +mn(C 2323 +C 2233 )+n£(C 33]2 +C 23n )+m£(C 22]3 +C 2312 ) 

m, 1, n are the direction cosines of propagation direction 

2.4 Phase and Group Velocity 


Phase velocity is defined as the velocity with which plane wave crests and troughs travel 
through a medium and is expressed as the ratio of the frequency of vibration and the wave 
number (i.e. the number of wavelengths per unit distance normal to the wavefronts). Group 



velocity, also known as ray or energy velocity, is defined as the velocity with which the 
energy of the wave propagates. 

An additional complication in anisotropic media is that the wave normal does not 
necessarily coincide with the propagation direction of the energy (ray direction). An 
arbitrary wave is dependent on time and spatial coordinates as / = f (cat — kx). In 
dispersive media, where the velocity is dependent on frequency, the acoustic energy of a 
narrow-band tone-burst propagates at the group velocity Vg. The phase velocity is simply 


and the group velocity is given as 



y _ d(D 


(2.24) 


(2.25) 
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Figure 2.2 Deviaiton between phase direction and ray direction 


In the symmetry directions only pure modes propagate and the phase and ray directions are 
necessarily parallel. This is because deviation from a symmetry direction in any direction 



must have the same effect as an equal deviation in the opposite direction. This also means 
that the velocities of pure waves are either maxima or minima for that particular mode. In 
any other directions the beam becomes skewed as illustrated in Figure 2.2. Needless to say 
that this skewing effect makes it more difficult to align the ultrasonic transducers for optimal 
inspection in highly anisotropic materials such as unidirectional composite laminates. The 
same effect is also responsible for the lateral shift of the transmitted ultrasonic beam through 
anisotropic plates that occurs even at normal incidence. 

2.5 Energy Flux and Group Velocity 

For acoustic waves propagating in elastic solids, the kinetic T and strain U energy densities 
are given as 


T = ~pu\ (2.26) 

U = ^C,„ m e u E, m , (2.27) 

Both T and U are scalar quantities. The total energy contained in a volume V is given by 

E = l(T + U)dV 

and the rate of change of energy with respect to time is given by 

£= dE, W + U) dv 
dt ly dt 

Using Eqs. (2.26) and (2.27) and the chain rule in the derivatives, Eq. (2.29) can be written 
as 


(2.28) 


(2.29) 
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a slight mathematical treatment and use of Gauss’ theorem will give 
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dV + ja ik udS k , 


(2.31) 


From equation of motion 


pii i 


da 


ik 


dx, 


= 0 , 


Eq. (2.31) reduces to 


(2.32) 


E = ja ik udS k 


Here by introducing the vector of energy flow density P, that is often called the acoustic 
Poynting vector, as follows 


P,=-aA- (2-33) 

P is directed towards the local direction of energy transfer. Finally, from Eq. (2.32) 

£ + |P.<iS = 0, (2.34) 

where dS denotes the elementary surface area vector that is everywhere normal to the 
surface S which encloses the volume V . 

Using the plane harmonic wave solution (2.9) the kinetic and strain energy densities can be 
given as 


T = ^pii 2 = ^ pco 2 A 0 2 sin 2 (k„x n - cot) 


(2.35) 


and 



(2.36) 


U=^C ijkl k J k,4p,p k sin \k n x n -cot) 


respectively, where it is exploited that the polarization vector is defined as a unit vector. 
From the equation of motion, 


pco 2 u — C.Jtk.u . 

* 1 ijkl j 1 m 


(2.37) 


Multiplying both sides of the equation by u, gives, 

pO) 2 u 2 = C iJkl k j k l uu i , (2.38) 

which leads to the conclusion that the kinetic and strain energy densities of a propagating 
plane wave are always equal 


T-U (2.39) 

Propagating waves exhibit spatial and temporal invariances that require that the kinetic and 
potential components maintain the same balance everywhere and all the time. In terms of 
energy density this means that the two components must be equal, i.e., they rise and fall 
together everywhere and all the time. It is a consequence of the Reciprocity Theorem, which 
requires that the transfer efficiency between two coupled forms in which energy can exist 
must be the same in both directions. Under these conditions, balance is feasible only if the 
kinetic and strain energies are equal. 

The time- averaged total energy density E and energy flow density vector P are given as 

~ 1 , 2 

E=-pcoA l (2.40) 

P' = ~2 Em Pj 


(2.41) 



Finally, the energy flow velocity Vg is defined as 




S' 


pco 


(2.42) 


2.6 Relation Between the Phase and Group Velocities 


The group velocity can be written in terms of phase velocity as 



pv 


(2.43) 


In order to fin d the relation between V g and V a vector Q is defined here as follows 


Q, = C ij, m n iP m Pj 


(2.44) 


group velocity can be written as 



(2.45) 


From the equation of motion (2.36) 

C IJlm n J n,p m p i -pV 2 =0. 


This can be written as 


Qn = pV 2 

Finally, from Eqs. (2.46) and (2.48), 

V n.=V 

& * 


(2.46) 


(2.47) 


(2.48) 


Or 



Vgcosy = y, 


(2.49) 


where y is the angle between V g and V. This relation indicates that V is the projection of V g 
in the direction of wave propagation 



Figure 2.3 Relation beween Phase velocity vector Group velocity vector 


2.7 Slowness Surface 

An extremely useful insight of the effects of anisotropy can be obtained by examination of 
the slowness surfaces of a material. The correct interpretation of these surfaces can aid in the 
prediction of beam profiles in anisotropic materials. 

The slowness vector of a wave, m, is defined as 



(2.50) 


Hence m has magnitude given by the reciprocal of the phase velocity. Since, for an 
anisotropic material, a wave’s phase velocity magnitude is dependent on the phase velocity 
direction, it can be said that |m| is a function of n, i.e. 


m=m(n) 


(2.51) 



If the value of m is plotted for all possible values of n, i.e. all phase velocity direction, then 
a three-dimensional surface mapped out by the end points of the slowness vectors, is called 
slowness surface. Since there are three modes of wave propagation in anisotropic materials 
there exists three slowness surfaces corresponding to the three wave modes. 

For an isotropic material each surface, corresponding to a wave polarization, will be a 
sphere since m is independent of n (i.e. phase velocity is independent of direction in 
isotropic materials) 

Figure 2.6 shows a slice through the slowness surface for transversely isotropic. There are 
several points to notice from the figure 

1. Each slowness surface is not a circle, indicating departure from isotropy. 

2. The tiree slowness surfaces, one for each mode, are of different shapes. This immediately 
indicates that no two wave modes will behave in the same manner with in the composite. 

A further use of slowness surface arises in depicting group velocity i.e. energy direction. 
Theory shows that for any phase velocity direction n, the normal to slowness surface for that 
particular value of n is in the direction of the group velocity. 


A 


> 



Figure 2.4 A general slowness diagram for an anisotropic material 


The usefulness of this unusual parameter than conventional velocity, lies in the nature of 
anisotropy itself. In many cases (e. g., refraction at a plane interface or diffraction from a 
periodic structure) the wave direction is determined by the wave speed For example, 
according to the well-known Snell's law, the refraction angles of all reflected and 
transmitted waves are determined by die incident angle 0/ and incident velocity V as 
follows 


sin<9 r _sin0 i 

V~ 


(2.52) 


where the subscript r represents any of the refracted waves. This simple explicit relation for 
6 r becomes an implicit one for anisotropic cases 


sin 6 r _ sin 0. 


(2.53) 


With slowness instead of velocity, the above condition becomes 

m r (0 r )sin 9 r = m t (6. )sin 6 i , (2.54) 

i.e., the projections of the slownesses of all waves are the same on the interface. This feature 
can be readily exploited to determine the refraction angles in a graphical way. At first the 
horizontal slowness projection m=m.(6 j )sm 9. is determined from the angle of incidence. 
Then, for each particular refracted wave the refraction angle is determined from the 
intersection point of the slowness diagram with the same vertical line representing m. A 
special case in which the vertical line does not intersect a particular mode’s slowness 
diagram then that mode will no longer exist in propagation 




Figure 2.5 representation of reflected/transmitted wave modes using slowness diagram 


2.8 Beam Divergence 


Beam divergence is measure of how fast the energy within the beam is spreading as the 
beam travels. For example, if wave has zero divergence then the spatial beam width remains 
constant along the beam path and the associated energy will be highly concentrated. 
Conversely, a large beam divergence leads to wide beams with a diffuse energy spread. A 
good measure of beam divergence is obtained by considering a small change in phase 
velocity direction, A0 P , and finding the associated change in group velocity direction, A6 g . 
The modulus of the ratio of these two gives beam divergence 


A6? 






(2.55) 


If D=1 (as is always the case for waves in isotropic materials), then beam widths are as 
expected for a given transducer. However, if D>1, then for a transducer whose frequency, 
dimensions, etc. are known, then angular speed of beam will be much larger then expected. 








This arises because adjacent layers in beam, of only slightly different phase velocity 
direction, have widely group velocity (or energy propagation) directions. Therefore, the 
beam spread from a conventional transducer when injecting ultrasound into anisotropic 
material will be much wider and hence more diffuse then would be case for isotropic 
material. Conversely, when EK1 it is find that an ultrasonic beam injected into anisotropic 
material will be narrower and hence more concentrated in energy, than it would be in an 
isotropic material. 


A good indication of beam divergence behavior as a function of propagation direction (i.e. 
phase velocity direction) can be obtained by examination of the slowness surface of a 
material. The local radius of curvature at any point on the slowness surface is a measure of 
beam divergence as this controls the rate of change of the surface normal with distance 
along the surface. Since the surface normal is the direction of group velocity, and distance 
along the surface is proportional to change in phase velocity direction (m), it gives 


n(0,)| = |m(e,)| 
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Here, r c (0 p ) is the radius of the curvature of the surface. 
The beam divergence is given by 


D 



(2.56) 


(2.57) 


A larger radius of curvature leads to small beam divergence and vice \ersa. Inspection of the 
slowness surface can give a good quantitative feel for beam divergence. 

To illustrate the beam divergence effects discussed above, let us consider file situation 
shown in fig 2.6, where a transducer fires ultrasound into a homogeneous transversely 
isotropic medium. By varying the angle of transducer beam the effect of probe direction on 
beam width can be illustrated. It is the shape of beam emerging at the bottom. To understand 



the location of the emergent beam, which is not always directly under the transmitter, 
consideration of beam skewing is required. 


Transducer 



Figure 2.6 Beam divergence and Beam skew 


2.9 Beam skewing 

The beam skew is defined as the difference between the group angle and phase angel. A 
positive skew angel signifies that the group angel is larger than the phase angel. Hence for 
all wave mode and phase angel combination having non-zero skew angel, in addition to 
change in beam profiles (due to beam divergence), the emerging beam’s positions will be 
altered from where it would be if the block were cf isotropic materiaLBeam skew,\|/, can be 
given as 

V = 0,-e r ( 2 - 58 ) 

giving 

dd p dO p 


(2.59) 





From this relation it can be clearly seen feat beam divergence is very low when rate of 
change of skew is ~ -1. Conversely divergence is high when the rate of change skew is large. 

2.10 Closure 

In this Chapter fundamentals of ultrasonic wave propagation in anisotropic media have been 
presented. The fundamental equation of wave propagation in inhomogeneous and anisotropic 
media called as CbristofFel’s equation has been discussed. Various terms associated with the 
elastic wave propagation in anisotropic media like phase and group velocity, slowness, 
energy flux and critical angle has been discussed. Beam divergence and Beam Skewing 
associated with elastic waves has been also discussed. 



Chapter 3 

Reflection and Refraction of Waves at a Planar 

Interface 


3.1 Introduction 

Nondestructive ultrasonic testing of inhomogeneous and composite materials is affected by 
several special features of wave propagation that arise from the strong anisotropy and 
inhomogeneity of these materials. The resulting complexity requires re-examination of old 
testing methodologies and development of new ones. One of the most fundamental 
phenomena in ultrasonic NDE is the reflection-refraction of ultrasonic waves at a plane 
interface. Even the simplest test procedure requires understanding of mode conversion and 
knowledge of elastic wave reflection and transmission coefficients and refraction angles. 
Reflection-refraction phenomena, while straightforward and well documented for isotropic 
materials, are much more complicated for anisotropic materials. In anisotrpic case Incident 
and reflected/refracted waves can no longer be thought of as purely longitudinal or 
transverse with the appropriate directionally constant wave speeds. In addition,- the direction 
of maximum energy flow does not, in general, coincide with the direction of the wave 
normal. 

3.2 Fundamentals of Reflection and Refraction 

The simplest situation is depicted in Figure 3.1a, where a wave encounters a boundary at 
right angle or normal incidence. The interaction only involves reflection of some of the 
wave and t ransmissi on of a portion, with the amount of energy in each part depending on the 
material characteristics. A more complicated situation may arise, particularly in solids, when 
the wave strikes at an angle, or at oblique incidence. What may occur, as shown in Fig. 3.1b, 
is that two types of waves are reflected for a single incident wave. This phenomenon is 



known as mode conversion, and is illustrated for the case of a longitudinal wave generating 
both longitudinal and shear waves. Yet another aspect is involved when waves encounter 
edges. Complex scattering and diffiuction of the waves may occur, similar to optics (figure 
3.1c). 

a) Incident Wave i « Reflection 






Figure 3. 1 Different type of acoustic wave interaction with material 
disconlmuties: (a) Normal Incidence (b) Oblique Incidence and Mode 
conversion (c) Diffraction and Scattering 


3.3 Reflection and Refraction in Isotropic materials 
3.3.1 Snell’s Law (Isotropic case) 

In the case of a plane interface between two isotropic elastic media in “welded" contact 
(rigid boundary conditions), imply continuity of all stresses and displacements across the 
interface. Figure 3.2 shows a schematic diagram of a plane wave with wavenumber kj 
incident on the interface at angle <j>,. The parallel lines with spacing equal to the incident 
wave length A.; correspond to equal-phase planes orthogonal to the incident plane. By 
definition, the wavenumber k* = 2 tc/A, is the magnitude of the wave vector k- The incident 
wave is converted at the interface into reflected and transmitted waves. The refraction angle 
of the transmitted wave is 4> r and its wave number is k- The wave fronts of the transmitted 
and reflected waves should match that of the incident wave along the entire interface. This 
follows from the requirement that the interface continuity conditions be satisfied not only at 
one particular point (e.g., at the origin), but everywhere along the interface. Any mismatch 
among the phases of the incident, reflected, and transmitted waves, that together must satisfy 
the boundary conditions, would lead to an inevitable variation of their relative contributions 
along the interface and thereby would prevent them from perfectly balancing each other over 
the whole boundary. 




Figure 3.2 A schemetic diagram of a plane wave incident on the interface 

From the universal phase matching requirement along the interface 

A sin <p. = A. and A sin (p r = A r (3.1) 

where A; and Ar are the wavelengths for the incident and transmitted waves and A, as shown 
in Fig. 4.1, is the spatial period of the wave motion along the interface. Equation (3.1) is 
Snell's law that can be also written in terms of the wavelengths 

a =_i_=_i_ > (3.2) 

sin 0,. sin (p r 

In terms of velocities 

V V 

V — — — — = — — — , (3.3) 

sin0j sin (p r 


In terms of wave number 


K~k i sin 0, = k r sin <p r 


(3.4) 



Here, Vi and V r are the phase velocities of the incident and refracted waves, respectively. 
Figure 3.3 shows the geometrical interpretation of Snell's law in terms of velocity vectors (a) 
and wave vectors (b), which can be summarized in the following way: 


a * hi 




Figure 3.3 Geometrical interpretation of Snell's law in terms of (a) Velocity (b) wave vectors 


1. The traces of the wavelengths (3.2) and velocities (3.3) of the incident, reflected and 
refracted waves on the interface are equal. The phase velocity trace V is shown in Fig. 
3.2(a). 

2. The projections of the wave vectors of the incident, reflected and refracted waves on the 
interface are equal. The wave vector projection K is shown in Fig. 3.2(b) (K = co/V). 

The same conclusion holds for multilayered media The projections on the interfaces of the 
transmitted and reflected wavenumbers always equals to that of the incident wave. In 
general at an interface, the incident wave is converted into longitudinal and transverse waves 
in such a way that Snell's law will hold for both types of waves, as shown in Fig. 3.4. Here, 
the incident wave kj produces four waves, namely a reflected longitudinal wave kn, a 
reflected transverse wave kn, a transmitted longitudinal wave k ]2 , and a transmitted shear 
wave kt 2 . 


According to Snell's law, all five (incident, two reflected, and two transmitted) waves have 
the same wave vector projection K on the interface. This is the basic concept on which 
reflection and refraction analysis of waves in the case of anisotropic and inhomogeneous 
materials can be based In anisotropic case the incident wave is generally converted at the 
interface into three possible modes of elastic waves and Snell's law is similarly satisfied 



Figure 3.4 Reflected and Transmitted wave vectors produced at 
an isotropic solid/solid interface by an obliquely incident wave 


3.3.2 Reflection/Refraction Analysis Using Slowness Surfaces 

The concept of the slowness surface has been introduced previously. Since in isotropic 
materials the wave velocity is independent of the propagation direction the slowness surface 
is a sphere with radius m = 1/V. As an example, Fig. 3.5(a) illustrates the case of a plane 
wave incident from a faster to a slower medium (mode conversion is avoided for the sake of 
simplicity), i.e., VI > V2. The wave vectors for the incident, reflected and transmitted 
waves are k;, ki, and k 2 , respectively, and <(>, is the incident angle. The wave vector of the 
incident wave originates on the slowness surface and is directed at angle <j), toward the origin 
as shown in the figure. Angles of the reflected 4> i and transmitted <(>2 waves can be found by 
equating the projections on the interface of all wave vectors involved 

Figure 3.5(b) illustrates the case of incidence at grazing angle Q; = 90°). The refraction angle 

V 

is clearly determined from the figure assin0 2 =-y-. The case of a wave arriving from a 

" T 



slow to a fast medium at oblique incidence is illustrated in Fig. 3.5(c). Here the retraction 
angle is larger than the incident one, 4*2 > <t>i. An important phenomenon occurs when (f), 
teach the critical angle <f) CT as shown in Fig. 4.4(d). At the critical angle the refraction angle 
<j >2 becomes 90° and the refracted wave becomes evanescent, i.e., an inhomogeneous wave 
that exponentially decays into the lower medium. In the absence of attenuation this wave has 
only reactive power and does not carry energy, however it may affect the phase of the 
reflected signal. 




c) f'i < f '2 


n < I- 2. $2 w 




Figure 3.5 Use of slowness diagram in finding the reflected and refracted waves 
in different cases 


There exist a “critical” incident angle which a wave makes with the traction free surface of 
an elastic half space and that as the incident angle passes through this critical value, the 
energy flux vector of the reflected' or refracted wave passes through grazing incidence and 
emerges as a surface wave. This wave is like a Rayleigh wave in that the amplitude of 
motion decay exponentially away from interface, particle motions have an elliptical orbit, 
and the wave travels with a constant phase velocity 



3.4 Reflection and Refraction at an Interface between two 
Anisotropic Media 

3.4.1 Snell’s Law (Anisotropic Case) 


Application of Snell’s law for anisotropic materials can be understand by considering a 
monochromatic plane wave incident on the interface between two anisotropic solids as 
shown in Fig. 3.6 the interface coincide with ?fX 2 plane and b is a unit vector parallel to the 

X3 


The displacement vector u, of an elastic wave can be given as 

u j = APj exp[t(V/ “ 6# )] (3.5) 

In terms of the previously defined slowness vector m, the displacement vector can be also 
written in the following alternative form 

u j = APj ex P ~ 0] (3.6) 



anisotropic media 


32 



Because of assumed rigid boundary conditions at the interface, it requires continuity of both 
displacements and tractions across the boundary. Then, continuity of the displacements and 
tractions at die boundary (x3 = 0) is given by 


and 



(3.7) 


a - 4 


3 

6 



(3.8) 


a=] a =4 


The incident wave is marked with superscript 0, and the three reflected and three transmitted 
waves are marked with superscripts 1; 2; 3 and 4; 5; 6, respectively. 


Using the relevant displacement- strain and constitutive relationships, the stress tensor can be 
expressed as 


= C ikjm ^ = ioJm m A oPj exp[iQ){m,x, - f)l 

777 

Substituting Eq. (3.6) into (3.7) the displacement continuity (3.7) can be written as 


(3.9) 


4)V,° exp (icamlx* )= p“ exp (icom“x“ )+ ^A^p* exp (ianfxf), (3.10) 

a-\ a=4 

Similarly Substituting Eq. (3.9) into (3.8) the stress continuity (3.7) can be written as 

m l c LjAp) exp (iaym'x,) 

= -'t<C , i3jm A^p a j exp (icomfx,) +tfn a m C" jm A 0 a p < ; exp (icom'x,). 


(3.11) 


where c'ugm and C n jkj ra , are the elastic coefficients of the first and second media. The 
boundary conditions must be satisfied not only for all times t, which is a direct result of the 
common term exp (iCQt) in all equations, but also at all points (xl, x2, x3 = 0) on the 

interface. Since the exponential functions exp(/Ct)7n“x / ) are linearly independent for a= 0, 
1, . . , 6 Eq. (3.10) and Eq. (3.11) will be satisfied only if all seven exponential functions are 



identically equal. Snell's law follows immediately from this simple feet since this condition 
requires that 


m)x l = m\x l = m]x l = m)x l - m l x l = m\x l — m\x l (3.12) 

Here the coordinate vector x lies on the interface plane therefore 


b-X = 0 (3.13) 

where b is normal to the interface as shown in Fig. 4.5. Another way of looking at Eq. (3.12) 
is that for any two different (3 and y values from the set (0,1, . • ,6) 

(m^ -m r )-x = 0. (3.14) 

Comparing Eqs. (3.13) and (3.14) shows that the difference m* 3 - m r is always parallel to b, 
Le., the plane defined by any two of the seven slowness vectors is always perpendicular to 
the interface. By definition, the slowness vector of the incident wave nf lies in the plane of 
incidence, therefore all the other surface vectors must also lie in the plane of incidence. 
Since the direction of the slowness vector coincides with that of the wave vector, it can be 
conclude that the wave vector of the incident wave as well as all the other wave vectors of 
the reflected and refracted waves lie in the same plane which is normal to the interface, i.e., 
the plane of incidence. Furthermore, 


-m r )xb = 0 or 

m^xb = m r xb, (3.15) 

i.e., all vector products of m° are equal, which is an alternative form of Snell's law. This 
relationship can be converted into a more familiar form by introducing polar angles <{) a 
between the interface normal b and the different wave normals n a for the incident wave and 
all refracted and reflected waves so that b. n a = cos (f> a . Then, according to the definition of 
the vector product 


m“ xb = m“sin <p a 


(3.16) 


O A 


Since b is a unit vector and oc=l ,2. . .,6 . 

Equation (3.16) expresses generalized Snell's law for arbitrarily anisotropic media. It is 
important to note that the slowness values m a in general depend on the angles 0 a due to 
anisotropy. Figure 3.7 illustrates the fundamental principle that slowness vector of the 
incident wave and all the other slowness vectors of the reflected and refracted waves lie in 
die same plane which is normal to the interface, i.e., the plane of incidence, and that they all 
have the same projection on the plane of the interface. The same conclusion holds for the 
projections of the wave vectors since k K = nfco. In terms of the wave vectors, we can 
summarize Snell's law as follows: 

> The wave vectors for the incident and reflected/refracted waves all lies in a single 
plane, "which also contains the interface normal and is called the incident plane. 

> All their projections on the interface are equal. 



Figure 3.7 Generalized Snell's law for arbitrarily 
Anisotropic media 


3.4.2 Reflection/Refraction Analysis Using Slowness Surfaces 

As concluded in the previous article incident and. all reflected and refracted waves lie in a 
single plane, i.e., in the incident plane. In a case where incident plane is xl x3 as shown in 



Fig. 3.8. In this “physical" coordinate system the slowness vector has only two 
components; and , while the components identically vanish. From Snell's law 
(3.16) the tangential components m a i in the plane of the interface are all equal to each other: 


m° = m\ 



(3.17) 


Here, rP i = sin<|) 0 , where n° is wave normal for incident wave, <|) 0 = <j>‘ is the incident angle, 
and V* is the phase velocity of the incident wave. Since n?i for the incident wave is known, 
the tangential components of the slowness vectors m“i (i.e., the projections of the slowness 
vectors on the interface) are also known. Equation (3.17)) can be re-written in terms of die 
corresponding wave vector components as follows 


k° =k]= = ft sin (f) r 


(3.18) 



Figure 3.8 Reflected and transmitted wave vectors produced at an anisotropic 
solid/solid interface by an obliquely incident wave 


From equation (4.19) only the tangential components m a i of the slowness vectors m are 
known and the normal slowness components m a 3 of the reflected and refracted waves are 
still unknown and must be determined. 



To find Hie slowness vectors for the reflected and refracted waves Snell's law and one 
additional equation is needed As shown in previously, the acoustic wave must satisfy the 
so-called Christoffel equation given as 

( C v u n j n t-p V2S ,Jp,= 0 - ( 3 - 19 > 

In terms of the slowness vector rr§ = n/V, this equation can be re-written as follows 

( C i jU m j m t -P 5 JP,= 0 (3.20) 

In order t> have a nontrivial solution for p : , the determinant of the coefficients this equation 
must vanish. 





By introducing a tensor T, as 


1 n=C lJk! mm k -pSi, 


(3.21) 


(3.22) 


Eq.(3.21) can be written as 


1^1 = 0 (3.23) 

In the chosen coordinate system, the plane of incidence coincides with the xl x3 plane, 
therefore me = 0 and Eq. (3.22) reduces to 

T a = [C ni/ (m,°) + (C,„ + C il3/ )m 3 m° + C^m\ ]- p8.„ (3.24) 

where mi 0 is used for all modes instead of mi due to Snell's law. This is a second-order 
polynomial in the unknown normal component nj of the slowness vector and can be written 
as 
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(3.25) 


1 \ l =a i i+b. l m 3 + d il m 2 m , 


where 


a „ = C nl ,(m°y -ps„, 
*„=(C m +C, 13 >,°, 

^ il ~ C'33 / ’ 


a u a h ’fyi b h ,d u d b 


(3.26) 


3.4.3 Geometrical Considerations on Reflection and Refraction 


As the projection of the group velocity vector V g on the wave direction n is equal to the 
phase velocity V 


V.n 

g' ‘ 


(3.27) 


The slowness vector m was defined as a vector parallel to the wave direction with a 
magnitude equal to the inverse of the phase velocity 


1 



i> 


(3.28) 


which leads to 


V, m = l. 


(3.29) 


Taking the exact differential of left side of Eq.(3.29) yields 

dV g • m + V, • dm = 0, (3.30) 


where 


/ 
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Here, the partial derivative of UK with respect to ft. can be calculated as 


dm 

7 

dn i 


dn, 


f n ^ 

A 



Hi.iL 

V 2 dn { 


(3.32) 


Substituting Eq.(3.32) into Eq.(3.3 1) we will have 


dm 8 . n. dV 

V p • dm = —V — ’-dn. =-V .(-£ )dn 

s gt *v_ i gi^ rr t/2 ' i 


dn 


V V 2 dn. 


1 rv dV ^ 
V K s ' dn } ' 


The right side of Eq. (3.33) is zero since from Eq. (3.27) dV /dn i = V gi , finally 


V -dm = 0. 

g 

This means that V g is always perpendicular to dm, i.e., to the slowness surface 

The slowness diagram is very useful in the determination of the reflection and refraction 
slowness vectors by graphical means. Let us assume that a quasi-longitudinal wave impinges 
at the interface between these two solids at an angle of incidence (j)j and draw a line in the 
incident wave direction through the origin of the slowness diagram of the first medium. The 
slowness vector m° of the incident wave is aiming at the origin from the intersection point of 
this line with the slowness curve of the quasi- longitudinal mode. According to Snell's law, 
the incident as well as all the reflected and transmitted waves will have the same tangential 
slowness vector component n?i, that is also indicated in Fig. 3.9. Now by drawing a vertical 
line on the opposite side of the origin at a distance n?i, the intersection of this line with the 
reflected and refracted slowness surfaces will determine the slowness vectors of the 
reflected and refracted waves. The group velocity vector (or ray direction) is normal to the 
slowness surface at these intersection points while the magnitude of the group velocity can 



be determined from the fact that the projection of the group velocity to the wave propagation 
direction is equal to the phase velocity, i.e., to the inverse of the slowness. 



t 



Figure 3.9 Analysis of wave 
reflection and refraction at an 
interface between two anisotropic 
materials based on the slowness 
surfaces 


3.5 Closure 


In this chapter refraction and reflection phenomena of ultrasonic waves at a planar interface 
has been discussed The use of slowness diagrams in refraction / reflection analysis has been 
discussed 
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Chapter 4 

Wave propagation in Transversely Isotropic media 

4.1 Introduction 

In this chapter the mathematical modeling of ultrasonic wave propagation in transversely 
isotropic media with arbitrary orientation has presented. For this solution of the equation of 
motion is derived for transversely isotropic media such as fiber composites or ideally fiber 
textured austenitic steels. The corresponding wave vectors characterize plane elastic waves, 
making especially possible a quantitative evaluation of the deviation of wave propagation 
direction and energy' flux, which is a basic characteristic of anisotropic materials. Reflection 
and refraction of plane waves at interfaces between two arbitrarily oriented transversely 
isotropic media is examined yielding an algorithm to trace the ultrasonic waves in anisotropic 
homogeneous or inhomogeneous material 

A unidirectional composite material may be modeled to be linear, elastic, homogeneous, and 
transversely isotropic with five independent elastic constant: Cun, C3333, C2323, C1212, and 
Ci 133. Figure 4.1 illustrates the symmetry characteristics of a transversely isotropic material. 
This symmetry class consists of an axis of symmetry normal to a plane of isotropy. 


jr> 



Figure 4. 1 Transversely isotropic unidirectional composite material 
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4.2 Equation of motion 


The plane wave equation in a linear elastic medium can be derived from the general equation 
of motion ( 2 . 7 ) and it can be expressed in tensorial notation as follows [ 19 ] 

(V • C • V) • u (R, co) + pco 2 u(R, cd) = -f (R, co), (4. 1 ) 

Where U is the displacement vector at any point, V is the gradient vector, p is the material 
density, f accounts for the volume force density, C is the elastic stiffness tensor and CO 

denotes the frequency. In the absence of body forces the equation can be modified as 

(V-C-V)-u(R,<) + p|^(u(R,0) = 0 ( 4 . 2 ) 

= Ot 

This equation of motion explicitly depends upon space R and time t 


The elastic stiffness tensor C for the general transversely isotropic materials has the following 
conventional relation. 
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(4.3) 


In terms of Lame’s elastic constant, elastic coefficients can be expressed as, Cim = 
A, + 2 fU L , C3333 =Aj. + 2 fl T , C2323 = / 4 r , Cni2 = P z . and C1133 = V . 

The stiffness tensor can also be written as 
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(4.4) 


C(a) = ^11 + ^ [(II) 1324 + (||) 1342 ] + (v - A T )[|aa + aa|] + 

(Ml - ^r)[(Iaa) 1324 + (aa|) 1324 + (|aa) 1342 + (aal) 1342 ] + 

[A, r + 2 fi T + X L + 2 ji L - 2(v + 2fi L )] aaa a, 

where unit vector a indicates the orientation of the material’s rotational sy mm etry axis, which 
will in the following be referred to as the “fiber direction.” I is the dyadic idemfactor, which 
can be decomposed, using dyadic products of the Cartesian unit vectors, according to 

J — 6 X + 6^6^ + 6 y Q y — 

accordingly, the tetrad I is the dyadic product of I with itself. The upper indidal notation 
indicates transposition of elements in the tetrads; 

The plane harmonic wave solutions of Eq. (4.1) with the stiffness tensor given by Eq. (4.2), 
i.e., the solutions to the homogeneous (f = 0) equation of motion, are in the form 


f 1 U U 

0 1 0 

0 0 1 


(4.5) 


u a (R,co) = Uu a ex V [JK a K • R], (4.6) 

where K is the propagation direction. 

In terms of slowness vector (m) Eq. (4.6) can be written as 

u a (R, co)-Uu a exp[ jcom a • RJ (4.7) 

where m a = ni/V a is the slowness vector. 

Solving the corresponding linear algebra problem leads to two dispersion relations 

ml SqL = p(B±(B 1 -4A)' /2 )/(2A), (4.9) 
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where 


A = ii L (X T +2ti T ) + (&-K) 2 (X L (l T +2fi r )-v(v 

+ 2fi L )) + (a •K) 4 [(/l i + 2}i L ){fi L -(Aj- + 2 fi T )) (4.10) 

+ (y + ^n 

B = h l +X t +2jj. t + (a-K) 2 (A z + 2jx L -(Ay + 2fi T )). (4.11) 

The abbreviations a=S, qS, qL denote Pure shear wave, quasi shear wave and quasi- 
longitudinal wave. 

The direction of energy transport or the group velocity can be given as 

(V f )„ = F„K + F„K„(rK)a-KK)-a (4.12) 


where 


F s = 


Bl-Bt 


B T + 0 i L -B T )(a-& 


K 


1 


qS>qL 


y a + 2(a-K) 2 z /4 ± 


F\ - (a -K) 2 F 2 
4b~ 2 - 4 a 


\ 


y A = A, (A, + 2\i t ) - v (v + 2 / 1 , ) 

z A = - (A r + /t r )) + (v + ^ £ )‘ 

^ = [pi - (K + 2 /0(A r + 2A£ r )][/i i - (A, + 2 - 
(v + fi L ) 2 (ji L + A r + 2fi T ) 


(4.13) 

(4.14) 

(4.15) 

(4.16) 


F 2 =F { +[nl -(A i +2/i i )(A r +2 J u 7 ,)][/i i - (A, + 2/iJ]- 

(v + + A £ + 2fi L ) 


A 

The group velocity can be given by the phase velocity V a = V a K plus an additional term, 
which depends on wave propagation direction K. and the fiber direction a. This shows tne 
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deviation of energy flow from propagation direction, which is known as beam skewing. It 

/V 

vanishes when K points in to symmetry directions. 

4.3 Refraction of plane waves at an Interface 

For understanding the refraction/reflection in the transversely isotropic medium the orientation 
of medium is taken such that the fiber direction is variable but parallel to the surface. Thus the 
unit vector a lies in the [a = (a x ,a y ,0)], as shown in fig.4.2. 



material 2 


material 1 


Y 


Figure 4.2 Fiber orientation , phase velocity and group velocity direction 
relative to cartesian coordinate system 


The plane wave solution (4.7) for the incident and the reflected/transmitted waves can be given 
as 


R (R,co) = U , u ex-pijm' R), (4.18) 

u *' ra (R,co) = U R ’ Ta u RJa exp(yW- ra • R), (4.19) 
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where for the problem in hand the slowness Vt \" Ta is to be determined. Assuming ideal rigid 

contact between two media requires the continuity of the particle displacement and the normal 
tractions according to 

uY+yiy*«u'“=xt/V\ (4.20) 

e-'l'+ie..-T' , “=Xe ! .T r “, (4.21) 

a a 

where T designates the stress tensor. 

Considering the wave propagation in the xz plane, Snell’s law requires all projection of the m 
vectors onto the interface to be identical 


1 Jta Ta 

m =m —m . 


and the y components of the slownesses to vanish 


m 


my - m y a = 0, 


(4.22) 

(4.23) 


with these conditions, the -components can be determined from the dispersion relations. 

Using dispersion relations (4.8) and (4.9) for refracted and reflected SH waves the m_ 
component can be given as 


= ( -b s ± -M -“&))/“* < 4 - 24 ) 

where 

a s =fl r + (fi L -fl T )al (4.25) 

b s =(l^ L -fi T )aa]s x (4.26) 

C s =[fl T +(V L -Hr)<K-P (4.27) 

For qS and qL waves, the m R ’ T<lP = S,L ) components of the slowness vector are defined as 

the roots of a fourth order polynomial, 
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(4.28) 


(mf) 4 (a q + b q a] + c q a z ) + (mf ) 3 • 2a x a z m x (b q +2 ca 2 z ) 

+ (mf ) 2 (m\ [2a q + b q (a x + a z ) + 6c q a x a] ] +d q + e q a z ) 

+ (mf ) • 2 a x a z m x (m 2 x (b q + 2 ca\) + e q ) 

+ m x (a q + b q a x + c q a x ) + m\ (d q + e q a x )+ f q = 0 

where 

a. q — fl L (X T + 2/l r ) 
b q =?i L (X T +2fi T )-v(v + 2fi L ) 

c g = (A + P-l )(Pi ~ (A + 2 Pt )) + (V + P L Y (4-29) 

d q =~p(P L + A + 2 Pr) 

e q =-p(A + 2 PL-(A + 2 Pr)) 

d q = p 2 


The Eq. (4.28) will give four values of m, values (two for reflection and two for refraction). In 
the present work HQR algorithm has been used for finding the roots of polynomial. The detail 
ofHQR algorithm is given in Appendix- A. 

To separate the roots corresponding to qS and qL wave the dispersion relation (4.9) is 
transformed with two conditions 


2 A 
P 



-B> 0 



-B< 0 


(4.3(f) 

(4.31) 


The corresponding group velocity vector for the refracted and reflected waves has to be point 
from the interface into the respective material. If a root’s pair is complex, a surface wave will 

occur and mf = 0 . 
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4.4 Closure 


In this Chapter the mathematical formulation ultrasonic wave propagation in transversely 
isotropic media with arbitrary orientation have been presented. Reflection and refraction 
analysis of a plane waves at interfaces between two arbitrarily oriented transversely isotropic 
media have been presented 
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Chapter 5 

Ray Tracing 


5.1 Introduction 

The growing need for the proper non-destructive testing of composites/composite laminates 
has raised considerable interest in studying elastic wave propagation in anisotropic media. 
The ray-tracing model developed in the present work deals with the propagation of ultrasonic 
energy through anisotropic and/or inhomogeneous materials. Energy is assumed to behave as 
bundles of rays, each ray acting independently from its neighbors. As mentioned in previous 
Chapters, fir anisotropic material, direction of energy propagation associated with each ray is 
generally not parallel to the phase velocity direction, the phase velocity defines the 
orientation of Hie phasefront associated with each ray. Furthermore, if a material is 
inhomogeneous then a ray path is likely to bend because of the continuous refraction as the 
ray experiences changing material properties. 

The model can be applied to the study of energy propagation in composite structures, to 
identify the regions that might be difficult to inspect because of poor penetration of energy in 
these regions. Such regions usually arise due to presence of inhomogeneities, or from the 
inherent nature of the material. As expected, the model has highlighted the complicated 
nature of energy flow in these structures. 

5.2 Assumptions 

1 . A plane wave solution of the wave equation has been assumed throughout. 
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2. It has been assumed that a geometrical construction for the ray path is valid i.e. no 
diffraction effects are allowed in passing through an inhomogeneous or homogeneous 
anisotropic medium. 

3. Boundaries between regions of different elastic properties are assumed to be perfectly 
smooth. 

4. Fiber composite materials are considered to be homogeneous anisotropic material. 

5. For a boundary between two media, one or both of which is inhomogeneous, it has been 
implicitly assumed that the stresses along the boundary can be calculated using the 
local elastic constants appropriate to the point at which ray impinges on the boundary. 

6. The waves are essentially monochromatic and continuous. Changes in ray parameters 

(i.e. phase and group velocities) due to adjacent rays in a beam have not been 

quantified. 

Validation of the Assumptions 

1. In general, for an ultrasonic beam not in the nearfield (N) of the probe, plane wave 

approximation is reasonable. For a probe crystal diameter d and ultrasonic wavelength 

A 


N-— (5-D 

42 

Using values appropriate to typical austenite inspection, d ~ 12.5 mm, X~ 1-5 mm, N ~ 
10-40 mm. 


For inspection of components whose dimensions are more then N, the assumption of 
plane wave propagation reasonable. 


2. For a geometrical construction to be valid, 


A« d 


5>-\Vc t - ^ '>>• is? gfcrsprai. 
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3jsflfccf 5?0 A 




(5.2) 
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where d is the dimension of any object which may cause diffraction. At lypical NDT 
frequency of a few MHz, the ultrasonic wavelength is around l-5mm Which is much 
smaller than the typical specimen or its defect dimensions. 

3. For material system of continuous parallel fiber, fiber spacing and the wavelength of 
propagating wave is of concern For the worst case it has been verified that the 
shortest wavelength is 350 pm (at 5 MHz) compared with the longest fiber spacing 
that can exist in a composite materials (=15 fim). Hence our material appeared to be a 
well-mixed continuum of anisotropic material with respect to practically used 
frequencies. So for the sake of simplicity we can leave the matrix - fiber interface in 
our analysis for reflection and refraction. 

4. For a boundary between two regions, neither of which has been formed by a weld 
process. Here surface roughness will be on the tens of micrometer scale and hence 
apparently smooth for a wavelength of few millimeters. 

5. In deriving the boundary conditions satisfied between two media, certain conditions on 
stress and displacement is specified. These were to be satisfied everywhere along the 
boundary. Here it has been implicitly assumed that the stress due to a wave of given 
displacement and slowness vector is the same at any point over the part of boundary 
covered by the beam width Strictly this is not so, due to change in elastic properties 

along the boundary. 

5.3 Two Point Ray tracing Method 

Two point ray tracing method consists of finding the path of that ray which will travel 
between two given points in a material. The two points generally represent transmitter and 
receiver position. For example, if a defect is situated in a inhomogeneous material as shown 
in fig. 5.1, then determination of the ray path that reaches the defect edge will require two 
point ray tracing. The rayl, emanating from transmitter in a straight path towards the defect 
will not reach the defect edge, because of refraction at interfaces and bending in 
inhomogeneous material. As shown in Fig.5.1, a 60° (angle measured from surface normal) 
ray aimed at the top of a defect is skewed, refracted and bent away from the required 
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endpoint So for determining the route by which, energy will reach the defect requires the 
solution of a non-linear boundary value problem, (the boundary values being the required 
initial and final ray position, the non-linearity arising from the distortion of ray path within 
the material). As there is no general analytical technique exists for solving this problem. 
Numerical iterative techniques have to be employed to find an adequate solution to this 
problem. Although the use of iterative techniques are computationally intensive, they are 
generally versatile in the sense that arbitrary inhomogeneity may be considered and as many 
arbitrary combination of different materials can be taken. 


Transmitter 



5.4 Shooting Method 

The shooting technique developed in the present work is a kind of ‘two point iterative ray tracing 
scheme’. Application of this shooting scheme allows calculation of the route that ultrasonic energ> 
takes from a source to the diffracting edge of a defect and to receiver. 

In ‘Two point iterative ray tracing” the start and end point (i.e. the transducer and the receiver 
locations) are specified by the user and the program determines the ray which passes through those 
two specified points. In two point iterative scheme there is no user intervention after the start and 
end points have been specified 
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In shooting technique developed here for a given transducer position, a conical shape fan beam 
of definite cone angle is sent in to the material and this beam is received at the other side. Here 
visualizing the path of the beam, the user can either increase or decrease the probe angle, such 
that rays can reach the desired receiver position. This shooting technique, coupled with the 
opportunity of user intervention, is far more efficient than its counterpart, namely, two point 
ray-tracing scheme, which requires large computation time because of its inherent iterative 
nature. Although the shooting technique also seems to be iterative, but is more efficient than 
the former as with this technique the user can interpret the effect of probe angle on ray path for 
a given anisotropic and inhomogeneous material. 

5.4.1 Applicability of the technique 

The success of the iterative technique relies upon the assumption that a sensible relation 
exists between initial slowness vector direction of a ray and its final location. This means that 
the change in initial direction of the slowness vector should lead to a corresponding change 
in final ray position, where the direction of these changes should not vary from one ray to a 
nearby ray. This will be true unless the ray paths cross each other in a nomuniform manner. 
Even if one spurious path is present, then it will slow down the convergence rate as well as 
generally impede successful solution. If large number of ray paths crosses each other in such 
a manner that the relationship between initial slowness vector and final ray position is not a 
smoothly varying function, then the technique will fail. 

The likelihood of misbehaving ray paths depends on the types of material used and on tire ray 
mode. The presence of many ray paths of this nature suggests that the energy propagation is 
incoherent and will be widely scattered. Furthermore, amplitude will be reduced because of 
the loss of coherence. The ray tracing approach is not appropriate in the vicinity of 
incoherent propagation. 

A cause of diffi culty, but not generally of failure, in the technique arises when small changes 
of slowness vector are accompanied by large changes in group velocity (energy direction). 
When this happens, the technique 'will be highly sensitive to the initial choice of slowness 
vector, as slight change in the initial direction of slowness vector will lead to large changes in 
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ray endpoint, reducing greatly the chance of reaching the desired endpoint within a 
reasonable number of iterations. 

5.5 Basic steps of the program 

The model is capable of predicting the propagation of ultrasonic waves within any 
component made up of one or more regions whose elastic properties and fiber orientations 
are known. The interface between two different regions is treated as a planar boundary such 
that reflection and refraction angles can be calculated using the elastic equivalent of Snell’s 
law. 

For tracing the path of different waves in a component, after specifying the size of 
component, the whole domain is discretized in to number of square elements. In the present 
work the size of element has been taken as 2 x 2 mm. In the case where an element has 
partially defect properties and partially base material properties, it is necessary to find the 
actual (‘edge’) interface normal direction for the sake of refraction and reflection analysis. 
For the sake of display purpose an element is shown to belong to one type material 
depending upon to which the major fraction of element area belongs as shown in Fig 5.2. 
Fig. 5.2 shows the meshing for a square block of side 200 mm. 




Defective elements 


Figure 5.2 Meshing of a square block of size 200mm *200 mm 



Ray paths axe calculated using a discrete stepping process. A ray travels a small distance in the 
current group velocity direction and the material properties at the start and end points of the 
step are found. The ray is then assumed to have crossed a boundary separating two regions 
with different material properties. This is implemented in the program by considering the 
refraction at a boundary between two regions. The elastic equivalent of Snell’s law is used to 
determine the new wavevector direction, which is dependent on the old direction and on the 
change in material properties across the interface (i.e. the properties at each end of a small ray 
step). The new ray parameters, such as group velocity and group angle etc, are determined 
using the theory of waves in anisotropic materials. 

Figure 5.2 is a flow diagram showing the basic steps involved in the program. In addition to 
checks shown in the flow diagram, there are numerous other error checks at each decision 
stage. There are however three principal steps in the algorithm: 

Step 1. Solution of the elastic wave equation for a general anisotropic medium leading to 
wave energy direction and associated particle displacement direction. 

Step 2. Stepping forward of the ray path in the direction of group velocity (energy direction) 
with, for inhomogeneous media, solution of appropriate boundary conditions across 
the step. 

Step 3. Calculation of refracted rays at a boundary between two regions of different material 
by applying appropriate boundary conditions. 

The combination of these principal steps, together with appropriate data storage and plotting 
routines, leads to simulation of ultrasonic wave paths. 

The program is written in Visual C++ having about 55 subroutines. In total there are about 
8000 lines of executable statements. A detailed documentation of the program and various 
subroutines is given in Appendix-C. 
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Figure 5.3 Flow diagram of Ray tracing program 


5.5.1 Case I: Homogeneous anisotropic media 

In the present model the medium is assumed to be non- dispersive. A ray traveling m some 
direction will continue in this direction until it reaches a boundary. Stepping forward is 
therefore relatively trivial as no new calculation of velocities is required after every step. It 
follows the ray through the medium, checking after each step that the ray is within the same 
region. This procedure continues until a change in region is detected. After a change in 
region has been detected, the exact location of crossing is found, together with the other 
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details such as direction of normal etc. This information, with the ray velocities and direction 
is then used for retraction analysis. 

/ 5.5.2 Case II: Inhomogeneous media 

For any inhomogeneous medium the elastic constants alters along a ray path. This leads to 
change in phase velocity and the group velocity. These effects will lead to a beam path, 
which is not a straight line but a curved one, the curvature of this path being dependent cm the 
variation of elastic constants over the region. 

In the program the whole region is divided in to number of elements and for each element 
elastic properties is specified. So for every new element in the path of the ray, the program 
compares the elastic properties of the new and the previous element and if there is any 
change in elastic properties, refraction at the interface between the new and the previous 
pixel will take place. Due to refraction there will be a alteration in wavevector direction and 
group velocity. The ray path is therefore composed of a set of lines showing alteration in 
direction. 

5.5 Closure 

In this Chapter various assumption used in the ray-tracing model and their validation have 
been presented. The basic steps involved in the algorithm have been also discussed. The 
difference in the analysis for homogeneous and inhomogeneous have been presented. 
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Chapter 6 

Results and Discussions 


In this Chapter, the results of wave propagation in inhomogeneous and anisotropic materials 

are presented. The results consist of ray paths developed by ray tracing model, slowness 

diagram, group velocity diagram for different modes and effect of the various fiber angle on 
various ray parameters. 

A two- dimensioned model has been developed which allows for the tracing of ultrasonic rays 

within a component, which may be inhomogeneous and anisotropic. The model also analyzes 

the reflection and refraction phenomena in anisotropic and inhomogeneous media. 

The effect of fiber angle and phase angle on wave propagation parameters like phase 
velocity, group velocity, group angle is presented. The effect of phase angle on beam 
divergence and beam skewing is also studied. The model also predicts the strength of signal 
received by the receiver of a given dimension. The results predicted by the model developed 
here are also compared with the experimentally obtained results. 

For the purpose of wave propagation analysis in anisotropic material, the following 
properties corresponding to Graphite/Epoxy composite are considered: . 

C] , = 145.8 GPa, C 33 = 13.5 GPa, C 44 = 3.4 GPa, Q 6 = 6.8 GPa, C I3 = 10.6 GPa and density 
=1.6 g cm' 3 

6.1 Wave propagation characteristic in Graphite/Epoxy 
Composite 

6.1.1 Phase and Group velocity 

In anisotropic materials the direction of wave vector (phase velocity) and the direction of 
energy flow (group velocity) are not the same and the magnitude of the phase velocity and 
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in the wave vector direction. So phase velocity is in general less then the group velocity. 

'Ihe variation o! group velocity / phase velocity with phase angle has been plotted for 
Graphitc/Epoxy composite with various fiber angles (Fig.6. 1-6.4). The behavior of the curve 
and the values ot group and phase velocity are quite matching with the experimentally 
obtained results by Degtyar et al [21] Fig. (6.5-6.6). 

Figure 6.1 shows the variation of group and phase velocity with the phase angle for quasi- 
longitudinal (QL) and quasi- shear (QS) wave in graphite epoxy specimen with 0° fiber angle. 

Figure 6.2 shows the variation of group and phase velocity with the phase angle for quasi- 
longitudinal and quasi- transverse wave in graphite epoxy specimen with 90° fiber angle. If 
can be very clearly seen that the behavior of 90° specimen is exactly reverse to that of 6 ] 
specimen which was expected. 

It is also clear that the velocity for quasi- longitudinal wave is generally largest when the 
wave vector (i.e. phase angle) is in the direction parallel to fiber direction. 

Figure 6.3 shows the variation of group and phase velocity with the phase angle for quasi- 
longitudinal in graphite epoxy specimen with 40° fiber angle. 

Figure 6.4 shows the variation of group and phase velocity with the phase angle for quasi- 
shear in graphite epoxy specimen with 40° fiber angle 

Here also it can be seen that for quasi- shear waves velocity is minimum when wave vector is 
parallel to the direction of fibers. 



Figure 6. 1 variation of Phase/Group velocity 
with Phase angle for Graphite/Epoxy composite 
having fiber angle 0 deg. 



Figure 6.3 Variation of QL Phase/Group velocity 
with Phase angle for GraphitefEpoxy composite 
having fiber angle 40 deg. 



Phase Angle (sfeg.) 

Figure 6.2 Variation of Phase/Group velocity with 
Phase angle for Graphite/Epoxy composite having 
fiber angle 90 deg. 



Phase Angte (deg.i 

Figure 6.4 Variation of QS phase/Group velocity 
with Phase angle for Graphite/Epoxy composite 
having fiber angle 40 deg 
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Figure 6.5 Experimentally obtained variation of phase velocity 
with phase angle for 90-deg. graphite epoxy composite C zO 



Figure 6.6 Expermentally obtained variation of group velocity with 

phase angle for 0 deg. Graphite/Epoxy Composite tan 
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Figure 6.7 shows the variation oi group velocity with the group angle for quasi- shear and 
quasi-longitudinal wave in graphite epoxy composite with 0° fiber angle. 



Figure 6.7 Variation of group velocity with 
group angle for graphite/epoxy composite 
having 0 deg. fiber angle 



Figure 6.8 Variation of group velocity with 
group angle for grapHie/epoxy composite having 
90 deg. fiber angle 



Figure 6. 9 Expermentally obtained variation of 
group velocity with group angle for 0 deg. 
Grp ahite/Ep oxy composite cnj 
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Group Velocity 


Fi c u!e .hews die \aiiadon of group velocity witH die group angle for quasi- longitudinal 
anti quasi- shear wave m graphite epoxy specimen with 90° fiber angle. 

1 i^uu. 6. ) shows the experimental obtained variation of group velocity with die group angle 
tor quasi- longitudinal and quasi- shear wave in graphite epoxy specimen with 0° fiber angle 

by Degtyar et al [21]. 

It can be seen that die behavior of 90° fiber is exactly reverse to that of fP fiber specimen. It 
can also be noted that the variation of quasi-shear waves of both 0° fiber and 90° fiber 
includes a cuspidal tegion. I his is due to fact that for waves at two different phase angle can 

have die same group velocity direction. 

Figuie 6.10 shows the variation of group velocity with die group angle for quasi- longitudinal 
wave in graphite epoxy specimen with 40° fiber angle. 

Here it can be seen that for tf to 90° phase angle variation the group angle varies ally from 
45.7 1 to 53.1° which shows that quasi- longitudinal waves in 40° fiber specimen has much less 
beam divergence than in 0° fiber or 90° fiber specimen. 



44 <8 50 52 54 


CrOLpAdQIS |®gt 

Figure 6.10 Variation of group velocity with 
group angle for graphite/epoxy composite 
having 40 deg. fiber angle 



Figure 611 Variation of QS group velocity with 
group angle for graphite/epoxy composite 
having 40 deg. fiber angle 
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Hune 6.11 shows the \ dilation of group velocity with the group angle for quasi- longitudinal 

"ave in graphite epoxy specimen with 40° fiber angle. 

6.1.2 Slowness Diagrams 

An extiemeh useful insight into the effects of anisotropy can be obtained by examination of 
die slowness sut faces of a material. Slowness diagram for graphite/epoxy composite 
specimen has been plotted. These diagrams are compared with the experimentally obtained 
by Spies [22). I he results obtained by the present model are quite similar to those obtained 

by Spies. 

Figure 6.12 and 6.13 shows the slowness diagram for quasi- longitudinal and quasi-shear 
waves in graphite/epoxy specimen with 0° fiber angle and 90° fiber angle specimen 

respectively. 



Figure 6.12 Slowness diagram for 0 deg. Figure 6.13 Slowness diagram for 90 deg. 

graphite/epoxy composite graphite/epoxy composite 
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Figure 6.14 and shows ihe slowness diagram for quasi- longitudinal and quasi-sHear waves h 

graphite epoxy specimen witii 40° fiber angle. 

It can be noted that gioup velocity direction is almost same between 0 - 90 ^ of pfiase angle. 



Figure 6. 14 Slowness diagram for 40 deg. 
graphite/epoxy composite 


By closely observing (Ire slowness diagram for quasi- longitudinal wave it can be seen that 
normal to the slowness diagram Is in the direction of wave vector for most of the time and as 
the direction of energy propagation (i.e. group angle) is given by the normal to the slowness 
surface. So from just observing the slowness surface of a material direction of energy 
propagation for a given wave vector direction can be predicted. 

6.1.3 Group velocity Diagrams 

Group velocity diagram for graphite epoxy composite specimen Has been plotted. These are 
compared with the experimentally obtained by Spies [22]. The results obtained by the present 
model are quite similar to those obtained by Spies. 
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Figure 6.15(a), (b) and 6.16(a), (b) shows group velocity diagram for quasi- longitudinal and 
quasi- shear waves respectively in graphite epoxy specimen with 0° fiber angle. 


m 



Figure 6.15 Group velocity diagram for graphite epoxy specimen 
with 0° fiber angle (a) QE wave (b) QS wave 




Figure 6.16 Group velocity diagram for graphite epoxy specimen 
with 90° fiber angle (a) QE wave (b) QS wave 
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6.1.4 Beam Skewing 


cm wan. piopagation and energy direction, known as beam skewing, is 

^ ^ ^ ' ltn< ^ ^ 01 Quasi- shear and quasi- longitudinal wave in graphite epoxy 

specimen w ith O' ' lilvr angle and 90° fiber angle. 



Figure 6. 17 Variation of skew with phase angle for graphite/epoxy composite 
having fiber angle 0 deg. 





CO 



Figure 6. 1 8 Variation of skew with phase angle for graphite/epoxy composite 
having fiber angle 90 deg. 
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Fioni the iinu obtained lot (1 fiber it can be predicted that between phase angles 16-30° the 

beam 4.cv. in.* is \er> high around 45-60°. 



Figure 6 1 9 Variation of skew with phase angle for graphite/epoxy composite 

having fiber angle 40 deg. 

Figure 6.19 shows variation of beam skewing for quasi- shear and quasi- longitudinal wave in 

graphite epoxy specimen with 40° fiber angle. 

It can be seen that the beam skewing vanishes in both the symmetry direction, i.e., along the 

fibers and normal to fibers. 


6.1.5 Beam Divergence 

Beam divergence is a measure of how fast the energy within beam is spreading as the beam 
travels. A plot between beam divergences against phase angle Is shown In fig. 6.20,6.21 and 
6.22 for 0, 90 and 40° graphite/epoxy composite. 
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Phase Angie (deg .) 

Figure 6.20 Variation of Beam divergence with phase angle for graphite/epoxy 
composite having fiber angle 0 deg 



Phase Angle (deg.) 


Figure 6.21 Variation of Beam divergence with phase angle for graphite/epoxy 
composite having fiber angle 0 deg. 







Figure 6.22 Variation of Beam divergence with phase angle for 
graphite/epoxy composite having fiber angle 40 deg. 


By observing the various plots beam divergence and phase angle it is seen that the beam 
divergence is highest when wave vector is in the direction parallel to the fiber direction. 

6.2 Ray propagation in Graphite/Epoxy Composite with defects 
6.2.1 Ray Paths 

For tracing the path of different waves in a component the whole domain (200 X 200 mm) is 
discretized into a number of square elements each of 2x2mm (Fig. 5.2). In the present 
analysis, three shapes of the defect are analyzed namely rectangular, circular and elliptical. 
M of these inputs are given in a visual interactive manner as shown in Appendix-B. 

The input parameter for tracing the ray paths are: 

> Base Material properties (i.e. elastic constants, density, fiber angle). 

> Defect location, size and shape. 

> Defect properties (i.e. elastic constants, density, fiber angle). 

> Initial Ray angle. 
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> Scanner position. 

> Receiver dimension and location 

> Wave mode to be traced. The option is provided to trace the same mode in the whole 
component or tracing different modes in different regions. 

Fig. 6.23 shows the basic layout of the software developed in the present work, which 
include component, ray paths, axes and various controls for analysis. 

The model can trace two types of beams, one is parallel beam and other one is fen rays 
(cone shape). For the present model the cone angle of 5° has been taken. The model can 
trace both the quasi-shear and quasi-longitudinal wave and also different modes in different 
regions. 



Fie idt Vtim Options Anafyts ttsfp 


d g? a i n> 


f 1 *.Bi 


■ 

& 


fS $ f :}.« j% 4® f? ft a S. 


JaTid: 



Material: Transverse^ isotropic media 
(fiber nxtgte 0 deg.) 

Defect: No Defect 
Ray An#e: 58 (teg, 

!§§ Qttasi-Lcm^idmal wave 

Mi Qaast-Shear wave 


Ready 


I pWHf 


Figure 6.23 Basic layout of the model 
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6.2.1.1 TteetasgMUar Defect 

Fig 6.24 shows both quasi-longitudinal and quasi-shear wave paths for 60° ray angle in 
graphite epoxy composite specimen having a rectangular defect. The size of the defect is 
40x 40 mm and the defect properties are taken to be Cn =60 GPa, C 33 = 30 GPa, C 44 = 3.7 
GPa, Cfi6 = 7.5 GPa, Ci 3 = 15 GPa and density =2.6 g / cm 3 



Figure 6.24 Ray paths in graphite/epoxy composite having a 
rectangular defect 

In this figure it can be that for a given phase angle 60° both quasi-longitudinal and quasi- 
shear wave have different energy propagation directions. Quasi-shear wave’s group angle 
37.30° while quasi-longitudinal wave’s is 44.47° 

6.2.1. 2 Circular Defect 

Fig 6.25 shows quasi-shear wave paths for 78° ray angle in graphite epoxy composite 
specimen having a circular defect of radius 40mm.The fiber angle of the base material is 0 
while the defect has 90° fiber angle. Scanner is positioned at 35mm from the origin (i.e. 
upper left comer ofblock). The defect properties are same as in 6.2. 1.1. 
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Figure 6.25 Ray paths in graphite/ejposy composite 
having a circular defect 


6.2. 1.3 Elliptical Defect 


Fig 6.26 shows quasi-shear wave paths for 65° ray angle in graphite epoxy composite 
specimen having a elliptical hole. The base material is 40° - graphite/epoxy composite. 


Hole’s minor radius =3 Omm, 

Hole’s major radius =50mm. 

Hole’s center = (110,130) 

Scanner is located at 20 mm from origin. 



Figure 6. 26 Ray paths in graphite/epoxy composite 
having a elliptical defect 
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^ beam divergence for a 40 ° fiber composite is very less than those having 
0°and 90° fibers. s 


6.2.1. 4 Mixed modes 

In this a ray is analyzed where in a quasi-longitudinal mode before ray is input, converted 
into quasi-shear mode in defect region and quasi-longitudinal mode after defect. Fig 6.27 
shows both mixed mode propagations for 67° ray angle in graphite epoxy specimen having a 
square defect. The side of the square defect is 50 mm. Scanner is located at 20mm from the 
left edge of block. 



Figure 6.27 Mixed modes wave propagation 

Figure 6.28 shows the beam type ray tracing in which a parallel beam of rays instead of a 
cone passes through the specimen. In this case base material has a rectangular defect of size 
80 x 40mm. 

It can be seen here that for 60° angle quasi-longitudinal wave passes through the defect but 
quasi shear waves just misses the defect 
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Figure 6. 28 Ray tracing in the form of a parallel beam 


6.2.2 Signal Strength at Receiver 

Signal strength is evaluated by means that for a given number of rays in a fen beam what 
percentage of rays send reaches a receiver position. Although there are various attenuation 
effects which came in to effect while a ray travels passes through a material but in the 
present work concentration is just made on the percentage of rays reaching to a particular 
position. 

Figure 6.29 (a) and (b) shows the signal strength at horizontal face for a given rectangular 
defect in a composite material. The receiver dimension is 25mm. Here it can be see that no 
rays can be collecting by putting the receiver on vertical face. 
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Figure 6.29 (a) Quasi-longitudinal ray path in a giaphite/epoxy composite having a 
rectangular defect (b) Signal strength at various receiver position 

Figure 6.30 (a) shows a path of a fan beam in a material having a circular defect. After 
refraction from the interface some of the rays goes to horizontal and some to vertical fece. 
Figure 6.30 (b) and (c) shows variation of the signal at horizontal and vertical fece 
respectively. 
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Figure 6.30(a) Quasi-longitudinal ray path in a graphite/epoxy composite having a circular 
defect (b) Signal strength at various receiver position on horizontal face (c) 
Signal strength at various receiver position on vertical face 

6.2.3 Critical Angle Phenomena 


As explained previously that there exists a “critical” incident angle, which a wave makes 
with the traction free surface of an elastic half space and that as the incident angle passes 
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through this critical value, the energy flux vector of the reflected or refracted wave passes 
through grazing incidence and emerges as a surface wave. 



Surface 

waves 


Figure 6.3 1 Critical angle phenomena 

Figure 6.31 shows 71° quasi- longitudinal ray paths in a composite having a rectangular 
defect. The rays which incident on the vertical face of the cbfect (because of the critical angle 
phenomena) converted in to surface waves and stops penetrating. 

6.9 Time of Flight 

The time taken by a ray to pass through a material is called time of flight While traveling the 
ray will get reflected from various interfaces It encounters in its path, lime of flight is a very 
useful data in the sense, which can be used for reconstructing the defect images by 
tomographic algorithm. Here in present work a fan beam of 50 quasi-longitudinal rays Is send 
in to the specimen and monitored on the other sides of specimen. 

Figure 6.32 (a) and (b) shows the 70° ray paths in a composite having a rectangular defect 
and the time of flight for each ray of the fan beam. The rays having incidence angle above 
70.8° do not reach the other face as they converted in to surface wave at the defect interface. 
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Figure 6.32 Time of flight data for a fan beam incident on a rectangular defect 

; igure 6.30 (a) and (b) shows the 60° ray paths in a composite having a circular defect and the time of 
light for each ray of the fan beam. It can be observed that between phase angle 61.8 to 62° there is no 
ay reaching to the opposite face and thus no time of flight recorded. 




Figure 6.33 Time of flight data for a fan beam incident on a circular defect 
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Chapter 7 

Conclusions and Scope for Future work 


7.1 Conclusions 

A two - dimensional numerical model has been developed which allows for the tracing of 
ultrasonic rays within a component, all or part of which may be inhomogeneous and 
anisotropic. From this model an understanding of physical reasons for the difficulties 
associated with ultrasonic inspection through composite materials has been gained. The 
model seems to be very useful to predict and understand the effects leading to beam 
distortion in anisotropic media. From the model developed in the present work following 
conclusions can be drawn. 

> The direction of energy flow in general does not coincide with wave vector direction 
when an elastic wave propagates in anisotropic media. Therefore, it is possible that 
the incident wave vector can be directed towards the interface but the energy can 
propagate outward from the interface. 

> Unlike isotropic case incident wave transformed into the three possible reflected and 
transmitted modes (i.e. pure shear, quasi- longitudinal and quasi shear). 

^ The wave propagation velocities (i.e. phase and group velocity) are maximum when 
wave propagates in a direction parallel to the fiber direction. 

> It has been found that if a specimen is scanned across a surface at constant beam 
angle, there is a possibility that ultrasonic ray energy may not penetrate certain 
regions. Therefore in such cases to ensure that the whole component can be 
examined, more than one probe angle may be needed. 

> The model is found to be a powerful tool in the numerical simulation of elastic wave 
propagation in transversely isotropic media. From the results it is found that the 
technique clearly models the special behavior of ultrasonic waves m such media. 
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7.2 Scope for Future work 


> In the present model although calculation of reflection angles is done but only the 
refracted waves has been plotted. The model can be extended so that it can also 
plot reflected waves. 

> Calculation of reflection and refraction coefficients can be incorporated • in the 
model. 

> Automatic change of step size for a ray path depending on then local rate of 
change of material properties. This should make the model more efficient when 
an inhomogeneous region has widely differing properties in different pails of the 

region. 
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Appendix - A 


This appendices will describe the numerical technique used for the solution of the forth order 
characteristic polynomial, the roots of which will give the z- component of slowness vectors. 
The polynomial is given by Eq. (4.28) 

( m f V( a q + b q a] + ca \ ) + (mf ) 3 • 2 aam x (b q + 2c a ] ) 

+ (mf )\m 2 x [2a q +b q (a 2 x + a 2 z ) + 6c q a 2 x a]] + d q + e q a 2 z ) 

+ (mf ) • 2aam x {m] (b q + 2c a] ) + e q ) (& 1} 

+ K (a q + b q a\ + c q a \ ) + m\ (d g +e q a 2 x )+f q =0 

The Eigen value method has been used for solving the above polynomial. 

Eigen value Methods 

The eigen values of a matrix A are the roots of the “characteristic polynomial” P(x) =det[A - 
xl]. This method is very efficient in terms of convergence and order of error. It can be easily 
verified that characteristic polynomial of the special mXn matrix 



is equivalent to the general polynomial 
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(a-3) 


m 

PW=|>/ 

i-O 

If the coefficients a t are real, rather than complex, then the eigen values of A can be found 
using the routines “balance” and “hqr”. It is a more robust technique, largely because of the 
fairly sophisticated convergence methods embodied in “hqr” routine. The method is to 
construct an upper Hessenbcrg matrix whose eigenvalues are the desired roots, and then use 
the routines balanc and hqr. 

Reduction to Hessenberg Form 

The strategy for finding the eigen system of a general matrix parallels that of the symmetric 
case. First the matrix is reduced to a simpler fonn, and then ail iterative procedure is 
perfomed on the simplified matrix. The simpler structure we use here is called Hessenberg 
form. An upper Hessenberg matrix has zeros everywhere below the diagonal except for the 
first sub diagonal row. For example, in the 6 X 6 case, the nonzero elements are: 

X X X X X X 

X X X X X X 

X X X X X 

X X X X 

XXX 
X X 

Such a structure can be achieved by a sequence of Householder transformations.Householder 
reduction to Hessenberg form is in fact an accepted technique. An alternative, however, rs a 
procedure analogous to Gaussian elimination with pivoting. This elimination procedure has 
been used in present work since it is about a factor of 2 more efficient than the Householder 

method. 
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Before the rth stage, the original matrix A has become A r , which is upper Hessenberg in its 

first r - 1 rows and columns. Hie rth stage then consists of the following sequence of 
operations: 

> Find the element of maximum magnitude in the rth column below the diagonal. If it 
is zero the stage is done. Otherwise, suppose the maximum element was in row r 

y Interchange rows r and r + 1. This is the pivoting procedure. To mak-R the 
permutation a similarity transformation, also interchange columns r and r + 1 . 

> For i =■ r+2, r+3, ,N, compute the multiplier 

«,>+■=— (a-4) 

a r+ >,r 

Subtract w, ir+] times row i+l from row i. To make the elimination a similarity 
transformation, also add n i r+1 times column i to column r+1. A total of N - 2 such stages are 
required. 

When the magnitudes of the matrix elements vary over many orders, try should be done to 
rearrange the matrix so that the largest elements are in the top left-hand comer. This reduces 
the roundoff' error, since the reduction proceeds from left to right. 

Balancing 

The sensitivity of eigen values to rounding errors during the execution of some algorithms 
can be reduced by the procedure of balancing. The errors in the eigen system found by a 
numerical procedure are generally proportional to the Euclidean norm of the matrix, that is, 
to the square root of the sum of the squares of the elements. The idea of balancing is to use 
similarity transformations to make corresponding rows and columns of the matrix have 
comparable norms, thus reducing the overall norm of the matrix while leaving the eigen 
values unchanged. A symmetric matrix is already balanced. 

Balancing is a procedure with of order F? operations. Thus, the time taken by the procedure 
of balancing, should never be more than a few percent of the total time required to find the 
eigen values. It is therefore recommended always to balance nonsymmetric matrices. It never 
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hurts, and it can substantially improve the accuracy of the dgen values c omp uted for a badly 
balanced matrix. 

The actual algorithm consists of a sequence of similarity transfo rm a tio ns by diagonal 
matrices D. The output is a matrix that is balanced in the norm given by summing the 
absolute magnitudes of the matrix elements. This is more efficient than using the Euclidean 
norm, and equally effective: A large reduction in one norm implies a large reduction in the 
other. 

The QR Algorithm for Real Hessenberg Matrices 

The QR algorithm with shifts can be given as 

Q,-(A,-*,I) = R. (a- 5) 

where Q is orthogonal and R is upper triangular, and 

a„,=r, q,+*,i 

= Q.-A,-QI (a-6) 

The QR transformation preserves the upper Hessenberg form of the original matrix A=A1, 
and the workload on such a matrix is 0(n 2 ) per iteration as opposed to 0(n 3 ) on a general 
matrix. As S — > 00 , A s converges to a form where the eigen values are either isolated on the 
diagonal or are eigen values of a 2 X 2 submatrix on the diagonal 

A key difference here is that a nonsymmetric real matrix can have complex eigenvalues. This 
means that good choices for the shifts k s may be complex, apparently necessitating complex 

arithmetic. 

Complex arithmetic can be avoided, however, by a clever trick. The trick depends on a 
lemma which states that if B is a nonsingular matrix such that 

B Q = Q H (a- 7 ) 
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where Q is orthogonal and H is upper Hessenberg, then Q and H are My determined by the 
first column of Q. 


The lemma is used in practice by taking two steps of the QR algorithm, either with two real 
shifts Ifc and kni, or with complex conjugate values k and k+i = kg*. This gives a real matrix 
Ash- 2 , where 


a„ 2 =q„, q, a, q: 


(a-8) 


The Q's are determined by 


A, -i,I = Ql R, 
a„,=q, a, q: 
A,.. -i,„I = QL R. 


j+1 1* n: j+ i 

Using Eq. (a- 1 0), Eq. (A- 11) can be rewritten 


Hence, if 


A„, - =Q» • Q„, • R„, Q 

M = (A,-i„,I)'(A,-M) 


Equations (a- 9) and (a- 12) gives 


R = Q M 


(a- 9) 

(a- 10) 
(a- 11) 

(a- 12) 

(a- 13) 

(a- 14) 


where 


Q=Q„i Q, 

R=R« R, 


(a-15) 
(a- 16) 


Hence Eq.(a-8) can be rewritten as 


A,<y=Q r a „ 2 

Thus suppose there is an upper Hessenberg matrix H such that 


(a- 17) 


A -Q T =Q T ■ H 


(a- 18) 
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where Q is orthogonal. If Q has the same first column as Q r (i.e.Q has the same first row 
as Q), then Q = Q and A t+2 = H . 

The first row of Q is found as follows. Equation (a- 14) shows that Q is the orthogonal matrix 
that triangularizes tire real matrix M. any real matrix can be triangularized by premultiplying 
it by a sequence of Householder matrices Pi (acting on the first column), P2 (acting on the 

second column) P n -i . Thus Q = P n _, P 2 • P, , and the first row of Q is the first row of 

Pi since P, is an (i — 1) X (z — 1) identity matrix in the top left-hand comer. Now Q must be 
find satisfying (a- 1 7) whose first row is that of Pi . 

The Householder matrix Pi is determined by the first column of M. Since A s is upper 
Hessenberg, Eq.(a-13) shows that the first column of M has the form >0,...,0], 

where 

P x = flu -a n (k s +k s+l ) + kk s+[ +a l2 a 2l 

( 7 i = a 2[ ( a u + a ii ~k s - & J+ i) (a-19) 

r \ ~ a 2i a 32 


Hence 

P, = l-2w, *w[ 


(a- 20) 


where wi has only its first 3 elements nonzero. The matrix P, ■ A s • P, is therefore upper 
Hessenberg with 3 extra elements: 


P, A • P, 


X 

X 

X 

X 


X X 
X x 
X X 
X X 


X X 
X x 
X x 
X X 
X x 
X 


X X 
X x 
X x 
X X 
X x 
X x 
X X 


(a-21) 
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This matrix can be restored to upper Hessenberg form without affecting the first row by a 
sequence of Householder similarity transformations. The first such Householder matrix, P 2 , 
acts on elements 2, 3, and 4 in the first column, annihilating elements 3 and 4. This produces 
a matrix of the same form as (a-21), with the 3 extra elements appearing one column over: 

x x x x x x 

x x x x x x 

x x x x x x 

x x x x x x 

X X X X X X 

XXX 
X X 

Proceeding in this way up to P n -i, it is seen that at each stage the Householder matrix P r has a 
vector w r that is nonzero only in elements r, r + 1, and r + 2.These elements are determined 
by the elements r, r + 1, and r + 2 in the (r - l)st column of the current matrix. Note that the 
preliminary matrix Pi has the same structure as P 2 ,. . . .,P n -i . 

The result is that 

P_, P 2 P, A, < < Pl.,=H (a-22) 

where H is upper Hessenberg. Thus 

Q=Q = P, P 2 P, (a- 23) 

and 

A, +2 = H (a-24) 


The shifts of origin at each stage are taken to be the eigenvalues of the 2x2 matrix in the 
bottom right- hand comer of the current A s . This gives 


k s "t” ^s+2 ^n-lyi-1 ^I,n 

^s^s+l ~ ® n,n ^ 


(a-25) 


Substituting Eq.(a-25) in Eq.(a-19) gives 
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Pi a 22{l( a nn Q U ~ A,*-l ]/« 2 l + a \l) 

g, =a 2I [a 22 -a u ~{a nn -a n )-(<W» -a,,)] (a-26) 

/■j — ^21^32 

here terms are judiciously grouped to reduce possible roundoff when there are small off- 
diagonal elements. Since only the ratios of elements are relevant for a Householder 
transformation, we can omit the factor a 2 1 from Eq.(a-26) 

In summary, to carry out a double QR step we construct the Householder matrices P r , r = 
l,....,n-l. For Pi we use pi, qi, and n given by Eq.(a-26). For the remaining matrices q, q, 
and r r are determined by the (r, r+1), (r+l,r+l), and (r+2,r- 1) elements of the current matrix. 

There are two possible ways of terminating the iteration for an eigen value. First, if a n>n .i 
become negligible then a nn is an eigen value. Then delete the n th row and column of the 
matrix and look for the next eigen value. Alternatively, q,-i , n -2 may become negligible. In this 
case the eigen values of the 2x2 matrix in the lower right-hand comer may be taken to be 
eigen values. Then delete the n th and (n - l) th rows and columns of the matrix and continue. 


93 



Appendix - B 


This appendix shows the various user interfaces provided for input the data like component 
size, material elastic properties, fiber angle, probe angle and wave modes etc. 

Basic View 

e Edit Vmm Options Analyil Help ^ __ __ _ 

□ & a m f ®i Vb *%! 


WE 




Component size 



Size of Specimen 


No. of Pixies in X- Direction ||" 


No. of PiKels in Y-Direction fo 


Note' Size of one pixel is 2 * 2 units 


OK 


Cancel ] 


Probe Angle 



Select Probe Afagel frost the List 




70 Degree 
fi 75 Degree 
p 60 Degree 
85 Degree 
SO Degree 
95 Degree 


OK | 
Cancel I 


Base Material 




Graphite Epoxy 1 

50.80 

Graphite Epoxy 2 

128.60 

Graphite Epoxy 3 

145.80 

MMC 

1 48.30 

CAS/SiC 

146.5 


25.50 


13.70 


.5 


165.50 


C66(G1 Fibre Angle 


5.8 


5.8 


-iSiisS: 

20 Degree 
25 Degree 
30 Degree 
35 Degree 
40 Degree 
45 Degree 
50 Degree 
55 Dearee 


Defect Shape 


jpefect Shape . ' ^ 
























dy the center of circular defect 

“ " V 


w 


ius jO 



osider Upper-Left corner as the Origin 

iaResiau? dimensioning accordingly 
inel Size 5*5. 


\ 


Location and Size 

Specify the Center of Elliptical 

x V Jo- 

Major Radius fo 

Minor Radius Jo” 


OK 


Cancel 



.Xj 


(land See— , 

!y the Upper-Left Corner | 

15 Y 15 j 


F 


OK | Carted' I 

ider Upper-Left corner as the Origin 
;eycur dimensioning accordingly. 


M 

I Pixel Nos. - — 

; Fran Pixel No. S' ToPixelNo. f5~ I 

'I I 3 



|~ Properties — : , - Parent Materia 


CltWN fo" ; 

C33fGPa) jo ““ j 

’ i 

C6S(GPa) W~ [ 

j If This is same as the 

■ . Parent Material Please 

1 check the Cheek boxJt 

j wiii Help u in visualizing 

j the Defects 

\ 

C13(SPa} JO 

; f” Parent Materiel 

i 

' i 

OWfGPa} [5 | 

. - ’ | 



~ rw uvecuun — — 

DertsSjr |0 j 

i 

Fiber angle jb 

, i 




,« -'feck 


OK 


Cancel 


Mae 
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Analysis for ray paths 



T"rrr"* 

& 




ill 


Sialic - - — __ 

f Check This if you want to trace same mod© in the whole region 

I Check This if you want region wise antyste (different modes in different regions} 


Cancel | 


Wave modes 


*j| 

Modes — : — * ~ i — — : j j 

I - ’ '' - ■ | f 

, Before Defect |rt Defect After Defect ] j 

j jQlogitudinal »| jo-Shear ' V| pffiSSBfil *| j j 


[ 0-Logiiudiri.j! ' 



lear - 


ok , : 


■ Cancel j 
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Scanner position 



XI 


Scanner Position — - 

> Scanner Position (X) jo 1 


Note: Assume Upper Left Comet of 
the Component as Origin 


° K I 

Cancel 1 


Receiver Information 




Dimension - 

; Receiver Dimension j0~ 

r-Face-- 1 

i ! 

! C Receiver on Horizontal Face ! 

J j 

; C Receiver on Vertical Face j 

Static — r~ ~ : — r_ j 

j f~ Points Marked ! 


ill 


-~0K~~| 

Cancel | 
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Appendix-C 


This appendix will describe the file structure of the program developed in the present work 
and it will discuss the various subroutines used in the program. 


File Structure 

In whole about 30 source code files and 25 header file are there in the program. Each file 
performs a specific task. The basic structure of the whole program based on these files is iven 
in die flow chart below. 
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Function of various files 

PixeLProp.cpp: It gathers the component size information from user and sends it for 
analysis 

MateriaLcpp: It gathers the base material propaerties from user and sends it for analysis. 
DefectjShape .cpp: It determines the shape ot the defect. 
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Defect_Shape_Circle.cpp: It gathers the size and location information of circular defect and 
sends it for analysis. 

Defect_Shape_Rect.cpp: It gathers the size and location information of elliptical defect and 
sends it for analysis. 

Defect_Shape_EHipse.cpp: It gathers the size and location information of rectangular defect 
and sends it for analysis. ^ 

Defect.cpp: It takes care if there is any particular pixel or group of pixels defective. 

Defect jShapeJProp.cpp: It collects the defect properties for any shape or pixel. 

Probexpp: It gathers the probe angle information from user and sends it for analysis. 

AnaIysis_Type.cpp: It finds out that user is interested to do fanbeam or parallel beam 
analysis. 

Mode.cpp: It gathers the wave mode for analysis. 

Ray.cpp: It stores the starting point etc. data for each ray. 

Scanner_Pos.cpp: It gathers the scanner position for analysis. 

Receiver.cpp: It collects the receiver size location etc. information. 

FirstJDoc.cpp: This file contains all the analysis subroutines. This is the mind of tire whole 
program. 

First_View.cpp: This file contains all the plotting subroutines. This is the heart of tire whole 
program. 

Receiver_curve: It calculates the data for signal strength at different faces and plots the 
curve. 
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Subroutines 

There are around 55 subroutines in the whole code. The main subroutines and their functions 
are as follows: 


mainjbeam: It analyses the input data for parallel beam of rays. 
main__fanrays: It analyses tire input data for fan beam of rays. 

christoffel_soln: It solves the Christoffel’s equation. It finds the group velocity and group 
angles for every step. 

ray_tracing: It calculates the ray path for foe given group velocity and angle. 

ray_tip _in_pixel : It calculates the pixel nos. in which tip and tail of pixel lies for a particular 
step. 

find_slowness_paralleI: It calculates foe parallel component of slowness vector. 

zrhqr, balance and hqr: This are numerical techniques used for solving foe fourth order 
polynomial. 

Iook_for_realroots: It calculates the normal component of slowness vector 
group_vel_spies: It calculates foe group velocity, 
stepjforward: It forward the ray by foe stepsize defined previously. 
check_for_crossing: It checks whether the ray has cross any boundary or not. 
red uced_stepsize: It sends back the ray to interface if ray has crossed any interface, 
output: It stores the coordinates of complete path of the ray and send it for plotting. 

OnDraw: It plots the ray paths. 
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